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13.  ABSTRACT  (Miximum  200  wordt) 

This  report  describes  a  modeled  system  comprised  of  an  infinite  shell 
with  both  inner  and  outer  fluids  subject  to  forced  harmonic  vibration  in  the 
longitudinal  (x)  and  circumferential  (0)  directions.  Two  excitations  are 
applied  to  the  outer  shell  surface:  normal  pressure,  P0,  and  longitudinal 
shear  stress,  Px.  The  equations  of  motion  are  derived,  from  first  principles, 
for  an  elastic  shell,  including  a  full  treatment  of  bending  and  rotatory 
inertia.  The  fluids  are  treated  as  linearly  elastic,  inviscid  media  and  are 
subject  to  the  wave  equation  in  three  dimensions.  The  stress- strain  relation¬ 
ship  for  the  shell  is  formulated  to  enable  the  treatment  of  isotropic, 
specially  orthotropic,  and  generally  orthotropic  composite  materials 
undergoing  plane  stress. 

Two  solutions  are  presented  for  the  equations  of  motion.  The  first  is  a 
harmonic  standing  wave  in  0  and  a  traveling  wave  in  x,  and  the  second  is  a 


14.  SUBJECT  TERMS 

Anisotropic 
Closed  Form  Solution 
Cylindrical  Shell 


Harmonic 

Mode  of  Propagation 
Orthotropic 


17.  SECURITY  CLASSIFICATION 

18.  SECURITY  CLASSIFICATION 

OF  REPORT 

OF  THIS  PAGE 

UNCLASSIFIED 

UNCLASSIFIED 

Structural  Damping 
Wave  Motion 
Wavenumber 


SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


15.  NUMBER  OF  PAGES 

144 


It.  PRICE  CODE 


NSN  7540-01-280-5500 


Standard  Form  298  (Rev  2-89) 
a,  *nv  via  f»-i» 

mioi 


UNCLASSIFIED _ 

security  classification 
OF  THIS  RACE 

13.  ABSTRACT  (Cont'd) 

harmonic  traveling  wave  in  6  and  x.  Modes  of  propagation  are  designated 
by  mode  number,  n,  which  defines  the  number  of  wavelengths  around  the 
circumference  of  the  shell.  Each  mode  of  propagation  consists  of 
components  from  the  radial,  longitudinal,  and  circumferential  (w,u,v) 
degrees  of  freedom  in  the  shell.  Although  the  system  is  capable  of 
supporting  an  infinite  number  of  modes  of  propagation,  only  the  first 
four  (beginning  with  n  =  0)  are  analyzed  here. 

A  three-dimensional  rendering  of  the  mode  shapes  provides  an 
understanding  of  the  deformed  shape  of  the  shell.  The  closed  form 
solution  for  the  magnitude  of  the  response  is  then  evaluated  for  three 
specific  shells  in  the  wavenumber- frequency  plane.  Results  are  presented 
in  the  form  of  transfer  surfaces  (i.e.,  magnitude  as  a  function  of 
wavenumber  and  frequency  for  a  particular  mode  number)  for  the  following 
cases: 

•  Interior  fluid  pressure  at  radial  position  r  normalized  by 
excitation  pressure  PQ, 

•  Strain  in  the  shell  wall  at  position  z  in  the  wall  normalized  by 
excitation  pressure  PQ,  and 

•  Interior  fluid  pressure  at  radial  position  r  normalized  by 
excitation  stress  Px. 

These  transfer  surfaces,  displayed  as  grey  scale  images,  illustrate  the 
dispersion,  as  well  as  cutoff  frequency,  group  velocity,  and  phase 
velocity,  of  the  branches  for  each  mode  of  propagation. 

The  problem  of  the  compliant  core  that  partially  fills  the  cross 
section  is  solved  for  the  n  =  0  mode,  and  these  results  are  also 
presented  in  the  form  of  transfer  surfaces. 
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NOMENCLATURE 


A 

Shell 

Total  cross-sectional  area  of  the  shell 

T 

=  Static  tension  per  unit  length  of 

A m 

= 

Cross-sectional  area  of  the  matrix  material 

circumference 

Af 

= 

Cross-sectional  area  of  the  reinforcing  fiber 

w,v,w 

=  Displacement  of  a  point  in  the  middle 

a 

= 

Mean  radius 

surface  of  the  shell  in  the  axial,  tangential, 

= 

Inner  radius 

and  radial  directions  (field  quantity) 

ao 

= 

Outer  radius 

u.v.w 

=  Components  of  the  amplitude  of  shell 

c; 

= 

Extensional  wave  speed 

motion  in  the  axial,  tangential,  and  radial 

r 

= 

Radial  coordinate 

directions 

X 

0 

: 

Axial  coordinate 

Angular  coordinate 

Fluid 

E 

Em 

— 

Young’s  modulus  for  an  isotropic  shell 

Young’s  modulus  of  the  matrix  material 

P  = 

Bulk  modulus 

Ef 

= 

Young’s  modulus  of  the  reinforcing  fiber 

cf  = 

ci  ~ 

Speed  of  sound  in  a  fluid 

Ey 

= 

Composite  Young’s  modulus,  1  direction 

Speed  of  sound  in  the  inner  fluid 

El 

= 

Composite  Young’s  modulus,  2  direction 

Ca  - 

Speed  of  sound  in  the  outer  fluid 

c 

= 

Structural  loss  factor 

= 

Particle  velocity  in  the  outer  fluid  due 

Cm 

:= 

Structural  loss  factor  for  the  matrix  material 

to  the  structure-induced  pressure  field 

C/ 

= 

Structural  loss  factor  for  the  reinforcing  fiber 

General  fluid  particle  velocity 

4> 

= 

Reinforcement  wrap  angle  relative  to  the 

= 

Inner  fluid  particle  velocity 

V 

= 

longitudinal  JC-axis 

Angular  coordinate  for  strain  or  displacement 

Pf  = 
Pj  = 

General  fluid  density 

Inner  fluid  density 

9 

. 

calculated  in  the  shell  wall  relative  to  the 
longitudinal  x-axis  of  the  shell 

Angle  formed  by  the  curvature  of  the  shell  in 

P  0  = 
K  = 

9/  = 

Outer  fluid  density 

Radial  wavenumber 

General  fluid  velocity  potential 

G12 

= 

the  r,  x  plane 

Composite  shear  modulus 

9,  = 
9i  = 

Inner  fluid  velocity  potential 

Outer  fluid  velocity  potential 

V 

= 

Isotropic  Poisson  ratio 

Po  = 

Normal  pressure  excitation  field 

= 

Poisson  ratio  of  the  matrix  material 

Px  ~ 

Axial  shear  stress  excitation  field 

v/ 

= 

Poisson  ratio  of  the  reinforcing  fiber 

Pe  = 

Circumferential  shear  stress  excitation  field 

Vl2 

= 

Poisson  ratio  (strain  in  the  2  direction  when 

Po  = 

Amplitude  of  the  normal  pressure  field 

V2I 

stressed  in  the  1  direction) 

Poisson  ratio  (strain  in  the  1  direction  when 

Px  = 
Pi  = 

Amplitude  of  the  axial  shear  stress  field 
Internal  fluid  pressure  field 

A 

— 

stressed  in  the  2  direction) 

Shell  thickness 

Ps  = 

V  = 

External  fluid  pressure  field 

Lame  constant 

1 'x 

— 

Mass  moment  of  inertia  about  the  x-axis 

Pf  = 

Lame  constant  (shear  modulus) 

I<b 

— 

Mass  moment  of  inertia  about  the  0-axis 

Jn  = 

Bessel  function  of  the  first  kind 

K 

— 

Fiber  misalignment  factor 

y„  = 

Bessel  function  of  the  first  kind 

p 

- 

Shell  density 

/„  = 

Modified  Bessel  function 

m 

n 

— 

Mass  of  the  shell  element 

Number  of  wavelengths  around  the 

= 

Hankel  function  of  the  first  kind 

Nij 

— 

circumference 

Shell  forces 

a  a 

Modified  Hankel  function 

Pressure  field  for  a  general  fluid 

°v 

= 

Stress  components 

xu  = 

Stress  tensor  for  the  fluid 

4 

= 

Strain  components 

Ekk  = 

Strain  tensor  for  the  fluid 

e/> 

= 

( i  &  j )  Represents  tensor  shear  strain 

8,  = 

Kronecker  delta 

Y/. 

= 

( l  ^  j)  Represents  engineering  shear  strain 

D,  E,  G,  H,  M  =  Fluid  velocity  potential 

neld  constants 

NOMENCLATURE 


Shell  &  Fluid 

c  =  Wave  speed 

cb  =  Breathing  wave  speed 

CO  =  Angular  frequency 

t  =  Time 

k  =  Axial  wavenumber  (2n/L) 

L  =  Axial  wavelength 

i  =  J— T 

Tr  =  Transfer  function  with  respect  to  radial 

excitation 

Tx  =  Transfer  function  with  respect  to  axial 
excitation 


Core 

ar  =  Stress  in  the  radial  direction 

Oq  =  Stress  in  the  circumferential  direction 

er  =  Strain  in  the  radial  direction 

Ee  =  Strain  in  the  circumferential  direction 

vf  =  Poisson  ratio 

Xc  =  Lame  constant 

Ec  =  Young’s  modulus 

uc  =  Radial  displacement 

uc  =  Radial  velocity 

rr  =  Core  radius 
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FORCED  HARMONIC  VIBRATION  OF  THE 
GENERALLY  ORTHOTROPIC  CYLINDRICAL  SHELL  WITH  INNER 

AND  OUTER  FLUID  LOADING 

INTRODUCTION 

The  analysis  of  wave  propagation  in  the  system  comprised  of  an  infinite  elastic  cylindrical 
shell  with  both  inner  and  outer  fluids  has  a  history  dating  back  to  the  beginning  of  the  19th 
century.  Thomas  Young1  (after  whom  Young’s  modulus  is  named)  derived  the  first  expression  for 
the  phase  velocity  of  long  wavelength  pressure  pulses  propagating  in  the  fluid  contained  within 
the  elastic  shell.  Korteweg  (as  described  by  Skalak2)  and  Lamb3  followed  with  similar  analyses  as 
efforts  were  made  to  understand  the  circulatory  system  of  the  human  body  and,  in  particular, 
pressure  pulse  propagation  in  blood  flow. 

The  phase  velocity,  cb  (given  in  this  report  by  equation  (72))  is  often  referred  to  as  the  Moens- 
Korteweg  velocity.  This  early  work  was  concerned  with  what  we  term  the  m>  component  of  the 
circumferential  ( n  =  0)  mode  of  propagation.  These  first  analyses  considered  only  the 
axisymmetric  excitation  and  response  of  the  shell/fluid  system. 

Warburton4  derived  the  response  of  the  isotropic  shell  with  a  limited  consideration  of  bending 
for  axisymmetric  and  nonaxisymmetric  excitation.  He  solved  for  the  resonant  frequencies  and 
wavenumbers  for  varying  shell  thickness-to-radius  ratios  and  for  various  circumferential  mode 
numbers.  Additionally,  the  isotropic  shell/fluid  system  has  been  studied  by  Gazis5  and 
Greenspon.6 

In  this  work,  we  are  interested  in  the  axisymmetric  and  nonaxisymmetric  response  of  the  shell 
that  is  “generally  orthotropic,”  which  means  that  the  principal  material  directions  are  not  aligned 
to  the  shell  axes.  Within  the  plane  stress  case  considered  here,  the  shell  is  also  anisotropic.  Thus, 
in  this  study,  generally  orthotropic  and  anisotropic  descriptions  of  the  material  are  the  same.  The 
spirally  wrapped  composite  tubes  widely  used  in  commercial  products  are  examples  of  generally 
orthotropic  shells. 

The  shell  is  treated  within  the  assumptions  of  plane  stress,  where  the  stress  in  the  r  direction  is 
assumed  negligible.  A  full  treatment  of  bending,  transverse  shear,  and  rotatory  inertia  is  addressed 
in  the  shell  development.  In  the  formulation  of  the  stress  resultants,  the  varying  circumferential 
length  of  the  hoop  fibers  is  retained.  Fluid  loads  are  calculated  at  the  inner  and  outer  shell  radii 
rather  than  at  the  middle  surface.  Every  attempt  has  been  made  to  retain  all  the  physics  that  are 
possible  within  the  scope  of  shell  theory  to  keep  the  solutions  valid  for  as  large  a  value  of 
thickness-to-radius  ratio  as  possible. 

The  solutions  for  the  displacements  of  the  middle  surface  of  the  shell  allow  for  the  calculation 
of  the  pressure  fields  in  the  inner  and  outer  fluids  as  well  as  the  stresses  and  strains  in  the  shell 
wall. 
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THE  PHYSICAL  MODEL 


The  model  derived  for  the  dynamic  response  of  the  system  to  harmonic  excitation  in  both  time 
and  space  will  consist  of  an  elastic  cylinder  that  is  expressed  as  a  membrane  or  bending  shell  of 
infinite  length,  which  contains  an  inner  fluid  and  is  immersed  in  an  outer  fluid. 

Three  forcing  functions  are  applied  to  the  cylinder:  pn,  px,  and  p0.  These  forcing  functions 
exist  as  pressure  fields  in  the  outer  fluid.  As  the  surface  of  the  cylinder  is  set  into  motion  by  the 
forcing  functions,  pressure  fields  p(-  and  ps  are  generated  in  the  inner  and  outer  fluids,  respectively. 

In  the  formulation  of  the  structure-induced  pressure  fields  p,  and  ps,  the  fluids  will  be  modeled 
as  compressible,  inviscid  media  that  will  not  support  a  shear  stress.  Since  px  and  p0  can  only  be 
imparted  to  the  cylinder  by  a  real  fluid  with  nonzero  viscosity,  this  approach  will  define  them  as 
given  harmonic  functions  at  the  outset. 


h 


Figure  1.  Physical  Model  Diagram 

First,  we  desire  to  know  the  deformation  of  an  arbitrary  point  of  the  cylinder  having 
coordinates  x,  0  and  distance  z  from  the  middle  surface.  However,  this  solution  requires  an  answer 
to  the  three-dimensional  elasticity  problem.  The  elasticity  problem  can  become  the  simpler  shell 
problem  if  simple  kinematic  relationships  are  established  between  an  arbitrary  point  in  the  shell 
and  the  corresponding  point  on  the  middle  surface  with  the  same  x,  0  coordinates.  When 
progressing  from  full  elasticity  to  the  shell  equations,  we  are  simplifying  the  elastic  behavior  in 
the  radial  direction.  The  following  assumptions  are  used  in  deriving  the  kinematic  relationships:7 
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1.  All  points  lying  on  one  normal  to  the  middle  surface  before  deformation  will  remain  the 

same  after  deformation. 

2.  For  all  kinematic  relations,  the  distance  2  of  a  point  from  the  middle  surface  may  be 
considered  as  unaffected  by  the  deformation  of  the  shell. 

3.  The  stress,  arn  in  the  r  direction  may  be  considered  negligible  compared  with  the 
longitudinal  stress,  a^,  and  the  circumferential  stress,  a^. 

4.  All  displacements  are  small  (negligible  compared  with  the  radii  of  curvature  of  the  middle 

surface),  and  their  first  derivatives  (the  slopes)  are  small  compared  with  unity. 

With  the  assumptions  stated  above,  we  have  simplified  the  elasticity  problem  to  a 
consideration  of  the  eight  stress  resultants  shown  in  equation  (1). 
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BENDING  SHELL 

From  the  shell  element  depicted  in  figure  2,  the  stress  resultants  (both  the  force  and  moments) 
are  derived  in  equation  (1)  on  a  per-unit  length  basis. 

Bending  stresses  arise  when  the  normal  and  shear  stresses  are  not  distributed  uniformly  across 
the  shell  thickness,  h.  When  this  condition  occurs,  a  moment  results  about  the  middle  surface  of 
the  shell.  The  1  +  z/a  term  in  the  stress  resultants  reflects  the  fact  that  hoop  fibers  of  different  radii 
have  various  lengths.  Accounting  for  a  shell  cross  section  that  is  trapezoidal  in  shape  causes 
NxQ  *  NQx  and  MxQ  *  MQx  in  this  development.  The  stress  tensor  is  unchanged  and  is,  as  usual, 
symmetric. 


MEMBRANE  SHELL 

When  we  choose  to  neglect  bending  and  twisting  moments,  the  shell  is  referred  to  as  a 
membrane.  The  membrane  forces  at  a  point  in  the  shell  represent  a  state  of  plane  stress  in  a 
tangent  plane  to  the  middle  surface.8  Three  stress  resultants  remain  from  equation  (1):  Nqq, 

and  N^.  In  the  membrane  case,  we  find  =  NQx  due  to  the  requirement  that  the  stress  resultants 
lie  in  the  plane  of  the  middle  surface.  The  trapezoidal  shape  of  the  cross  section  is  neglected, 
which  means  that  the  z/a  term  is  dropped  from  equation  (1).  When  this  is  done,  it  is  clearly  seen 
that  Nqx  =  N&.  Equation  (2)  gives  the  relationship  between  the  actual  stress  and  the  for  the 
membrane  shell: 


(2) 
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STRESS-STRAIN  RELATIONS 


A  shell  without  reinforcement  is  by  definition  homogeneous  and  isotropic.  Such  a  material  has 
the  simplest  stress-strain  relationship  and  can  be  defined  by  two  material  property  constants. 
Often  an  axial  reinforcement  is  embedded  in  the  shell  to  provide  an  increase  in  strength.  When  the 
principal  material  directions  are  aligned  with  the  axes  of  the  structure,  the  material  is  termed 
“specially  orthotropic.”  A  specially  orthotropic  composite  obeys  the  second  simplest  constitutive 
law,  and  four  independent  material  constants  are  required.  In  the  specially  orthotropic  case,  the 
normal  stresses  do  not  depend  on  the  shear  strain  as  we  can  be  observe  by  inspecting  the  Q  matrix 
in  equation  (3)  and  noting  the  location  of  the  zeros.  If  the  principal  material  directions  are  not 
aligned  with  the  coordinate  axes  of  the  structure,  we  would  then  have  a  fully  populated  Q  matrix. 
Such  a  case  would  be  represented  by  a  spirally  wound  reinforcement.  This  situation,  where  the 
principal  material  property  axes  are  not  aligned  with  the  coordinate  axes  of  the  shell  (figure  3), 
results  in  the  “generally  orthotropic”  lamina. 

MEMBRANE  SHELL 

Specially  Orthotropic 

Equation  (3)  gives  the  desired  relation  between  stress  and  strain  for  the  specially  orthotropic 
case  of  plain  stress:9 
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Using  equation  (2)  and  noting  that  Cy  =  hQjj,  we  rewrite  equation  (3)  in  the  following  form: 
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Isotropic 


Equation  (3)  will  collapse  to  the  isotropic  result  when  we  set  v12  =  v2i  =  v  and  £,  =  £2  =  E: 
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A  connection  between  strain  and  displacement  is  required  within  the  shell.  Flugge10  provides  a 
detailed  development  of  the  strain  displacement  relations,  which  are 


du 


E*x  dx 


1  ,5v  . 

"  a(dO  h 


dv 

dx 


Vxe  ~  ' 


1  du 

add 


(6) 


Caution  should  be  exercised  when  we  work  with  the  tensor  representation  of  shear  strain.  There  is 
a  difference  between  engineering  shear  strain  and  tensor  shear  strain,  where  y,y  =  2e,y.  The  use  of 
one  over  the  other  will  change  the  terms  in  the  stiffness  matrix  Q.  The  values  of  the  indices  (i.e., 
11,  12,  16,  26,  and  66)  relate  to  the  position  in  the  full  6,6  tensor  representation. 

Generally  Orthotropic 

The  stress-strain  relations  for  the  generally  orthotropic  lamina  are  obtained  by  rotating  the 
stresses  from  the  specially  orthotropic  case  through  the  angle,  0.  The  resulting  relations  are  given 
by  equation  (7)11  as 


(7) 
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Figure  3.  Material  Property  Versus  Shell  Coordinate  System 
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In  equation  (7),  the  bar  over  the  indicates  transformed  stiffnesses,  given  below  as 

C 


0 11  =  £?nCOS4<j)  +  2(Ql2  +  2Qfa)  sin2<J)cos2<))  +  £)22sin4<f>, 

0 12  =  (<2ii +Q22~4^66)  sin2<))cos2<t)  +  0i2(sin4<t>  +  cos4<)>)» 

022  =  £?nSin4<|)  +  2(@12  +  2@66)  sin2 <J> cos2 $  +  Q22 cos4 <(>, 

016  =  (011  -012 -2066)  sin<t>c°s3<}> -H  (012 -022  +  2066)  sin3 <)> eos 0, 
026  =  (011 -012 -2066)  sin34>c°S<j)+  (012 -022  +  2066)  sin<l,cos3^ 
066  =  (011  +022~2012-2066)  sin2<t)cos2<()  +  Q^i sin4 <j>  +  cos4 <j>), 


From  equation  (2)  for  the  membrane  shell,  we  write  equation  (7)  in  the  following  form  noting  that 

Cij  =  hQij : 


Cll  C12  ^16 

C\ 2  C22  ^26 
^16  ^26  ^66 


(9) 
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BENDING  SHELL 
Generally  Orthotropic 

As  will  soon  be  seen,  deriving  the  stress  resultants  for  a  bending  shell  is  considerably  more 
involved  than  for  the  membrane  shell.  Beginning  with  the  assumptions  governing  the  kinematics 
of  the  shell  (stated  in  the  beginning  of  this  sectio’ ),  we  arrive  at  the  following  strain  displacement 
relations:12 
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We  use,  as  before,  equation  (7)  for  the  stress-strain  relationship  for  the  rotated  system: 
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Deriving  expressions  for  the  stress  resultants  in  terms  of  the  displacements  u,  v,  w,  and  their 
derivatives,  we  substitute  equation  (10)  into  (11)  and  then  into  equation  (1).  Now  the  integrations 
can  be  performed,  and  we  arrive  at  equation  (13),  where  the  stress  resultants  are  in  the  required 
form.  The  rigidities  are  grouped  together  in  equation  (12)  and  the  stress  resultants  are  listed  in 
equation  (13): 
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Appendix  A  provides  one  complete  example  of  stress-resultant  integration  along  with  a 
discussion  of  the  simplifying  assumptions  that  were  necessary  to  arrive  at  the  final  form  of 
equation  (13). 
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EQUATIONS  OF  MOTION 


MEMBRANE  SHELL 

Figure  1  is  the  physical  model  diagram  for  the  infinite  shell  of  thickness  h  and  midradius  a. 
The  displacements,  u,  v,  and  w  are  in  the  x,  0,  and  r  directions;  pQ  is  the  excitation  pressure;  and  ps 
and  pj  are  pressure  fields  that  result  from  the  motion  of  the  shell.  Consider  the  incremental  volume 
of  a  shell  that  is  dx  long  by  ad©  wide  with  thickness  h  (see  figure  4).  The  equations  of  motion  for 
the  system  are  determined  by  summing  forces  and  setting  the  result  equal  to  the  mass  times  the 
acceleration  in  the  three  principal  directions.  The  are  forces  per  unit  length.  The  pressure  loads 
p0,  ps,  and  pi  have  units  of  force  per  unit  area. 


Figure  4.  Membrane  Shell  Element  FBD 


From  first  principles,  the  equations  of  motion  are  derived  for  the  shell.  Summing  forces  in  the  x 
direction  gives 
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Figure  5  illustrates  the  membrane  shell  element  in  the  r,  0  plane.  When  we  sum  the  forces 
acting  in  the  radial  direction,  Nqq  and  the  static  tension,  T,  in  the  axial  direction  must  be 
considered.  In  the  undeformed  state,  the  shell  is  parallel  to  the  x-axis  with  no  component  of  T  in 
the  radial  direction.  However,  when  the  shell  undergoes  deformation,  as  shown  in  figure  6,  a 
component  of  force  resulting  from  the  static  tension  arises  in  the  radial  direction.  This  static 
tension,  T,  is  much  larger  than  the  forces  resulting  from  the  applied  acoustic  pressure  loads  and 
may  therefore  have  a  significant  effect  on  the  structure,  depending  on  its  relative  magnitude  as 
compared  with  the  other  forces  involved  in  the  radial  direction.  For  this  reason,  static  tension 
must  be  included  in  the  formulation. 


Nqq  must  be  separated  into  components  in  order  to  sum  forces  in  the  radial  direction.  Small 
angle  approximations  will  be  used  because  9  is  small;  therefore, 

cosO  =  1,  sin0  ~  0  • 


Figure  5.  Membrane  Shell  Element  FBD  in  the  r,  9  Plane 
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Figure  6.  Membrane  Shell  Element  FBD  in  the  r,  x  Plane 


Thus,  we  find,  that  NQQr  =  Aeed0/2  and  Tr  =  Td(p/2.  The  elemental  length  in  the  axial  x-direction  is 
equal  to  r{d<. p.  Summing  forces  in  the  radial  direction  while  referring  to  figures  5  and  6  yields 
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Assuming  that  the  slope  of  the  surface  is  small,  the  following  form  for  the  radius  of  curvature  is 

1  d2w  dNeed02rld<p 

used:  —  =  — z.  Neglecting  the  higher  order  term  ^  - =- — ,  we  have  the  following 

r,  2 

equation  of  motion  in  the  radial  direction: 


Summing  forces  in  the  0  direction,  we  have 


m  =  pahdxdQ, 


-Needx  + 


3/Vee 

(■ Nee  +  90  dQ)dx  +  pedxadQ  -  NxQadO  +  ( Nxe 


dN. 


xQ 


,a2v 


+  ^-  dx)adO  =  p  ahdxdO-pr, 
dx  K  dt2 


which  can  be  reduced  to 


la/vee  dNxO  ,d2v 

a3e  +K  = 


(16) 


Summarizing  the  equations  of  motion  in  the  three  directions,  we  obtain 


1^9, 

Mxx 

~\  2 
,  3  u 

adO 

+  dx 

+  Px 

=  (V 

1  a^ee 

,  32v 

a30 

+  dx 

+  Pq 

=  PV 

. 

rd2w 

,d2w 

~Ps~ 

—  +  j 

a 

r - r- 

dx2 

=  PV 

(17) 

(18) 

(19) 


Equations  (17)  through  (19)  give  relationships  between  the  applied  pressures  and  the 
corresponding  displacements  and  stress  resultants.  We  will  proceed  in  a  direction  that  gives  three 
equations,  where  the  only  unknowns  are  u,  v,  w,  and  their  derivatives.  By  means  of  equations  (17), 
(18),  (19),  (6),  and  (9),  we  can  now  write  three  second-order  linear  partial  differential  equations 
for  the  displacements  of  the  cylinder. 


The  first  step  substitutes  equation  (6)  into  (9): 


..  ^  3k  ?  1  ,3v  „  -  ?v  .  ldn. 

Nxx  -  Cli^;  +  C12-  (3fi  +  vv) 


N, 


ee 


p  du  p  1  .3v  ,  .  ,  p  ,dv  ,  1  du. 

-  c<2S]r+c22-(3S+«'i  +c26(3i+-3S), 


p  du  .  p  1  ,3v 


N*  =  C16^  +  C26^(^  +  h-)+C66(^:  + 


3v  .  1  du. 


l63rl'“5'56T’v;T''«'3rS5e) 


(20) 
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Differentiating  equation  (20)  where  needed  yields 


~  d2u  p  1 

.c 

>2v 

3h\ 

3* 

'V12* 

(303jc 

+  K> 

+  ^16 

1  | 

f32v 

3vt^ 

30 

C'23xM  +  C2i 

a 

l302 

+  30; 

+  ^26 

p  d  2U  ~ 

1 

r32v 

dw' 

30 

C|6a^e+C“ 

a 

,302 

+  30; 

+  C66 

p  3  2«  p  1 

.3 

l2v 

3w> 

3x 

163?+  26  * 

303* 

+s> 

+  C66 

a303x  j 

fd2v 

1  32w  ^ 

1^3x30 

0302/ 

fd2v 

1 32n  ^ 

[3*30 

+  *302/ 

(d2v 

1 

1  32m  ) 

— ^  + 
l.3x2 

cz303x  J 

(21) 


Substituting  equation  (21)  into  equations  (17),  (18),  and  (19)  gives 


-c  d  “  ,  r  1  ^3w^  ^  1^  (d2v  1  32a 

a  169*3e  26«2  Ue2  56  J  a  +  <*302J 


+  c  a2"  + 
l+c"a7+ 


C  1 A  ,  ^  ,  r  P2v  ia2M  >1 
Cl2a(aedjc  ajc}  Cl6Uv2  <*303* 


+Pjc  = 


,  32m 
*3? 


(22) 


Cl2«a]ca0+C22^[^2+ae 


+  C  1  p2v  ,  I^h)  o  d2u 

26 a  [3*30  aae2J  l6dx2 


p  1  ,d\  .  3h\  . 


a  303.x  dx 


)+c, 


66 


a2v  +  ia2M 


aaeax 


+/?e  = 


,a2v 

•v 


(23) 


n  „  r.  ^  ia«  r  1  \  ^  1  ,  l&u\  .  rd2*  ,  a 

/>,  P»-',--C'^S-C22?(5S  +  »')-C26-(^+-3g)+7-^  =  pA^-. 

The  following  harmonic  excitation  will  be  applied  to  the  shell  in  the  r  and  ^-directions 
p0  =  />ocosn0e'(^-®f),  px  =  P  cos/i0<?,(**-o),) . 


(24) 


(25) 


Normal  to  the  surface  is  the  pressure  field  pa  and,  in  the  longitudinal  direction,  the  shear  stress  px. 
Both  P0  and  Px  are,  in  general,  functions  of  wavenumber  and  frequency;  however,  they  have  been 
expressed  here  in  a  simpler  form  as  constants  since  all  the  results  presented  will  be  normalized  by 
P0  or  Px  and  only  the  structures  response  to  a  unit  pressure  excitation  is  considered.  It  should  also 
be  noted  that  P0  and  Px  are  amplitudes  of  the  components  of  each  mode  of  propagation  and  are  in 
general  different  for  each  mode  n. 
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(26) 


The  solutions  to  equations  (22),  (23),  and  (24)  are 
u  =  (/cos  («0)  v  =  Vsin  («0)  e  ^kx~a,\  w  =  W co&  (nB)  e‘ (kx 

BENDING  SHELL 

Figure  7  depicts  the  forces  acting  on  an  elemental  area  of  the  shell.  In  this  case,  we  include  the 
transverse  shear  forces  that  can  now  be  supported  by  the  bending  shell  model.  Proceeding  as 
before,  we  sum  forces  in  the  three  coordinate  directions  to  obtain  three  equations  of  motion  for  the 
middle  surface  of  the  shell. 


Figure  7.  Bending  Shell  Element  FBD  Forces  Only 


IS 


Starting  with  the  r  direction,  we  have 

♦t 


a2w 


aw 

F  =  m~, 
dt2 


m  =  pahdxdQ, 


dQ 


0r 


de 


dN, 


ee  dQ 


dQ)  dx 


(0Qr  Qer  0q  dQ)dx+  ( Qxr  Qxr  dx)adQ  NQQ  ^  dx  +  ^ 

d2w 

+  pradQdx  =  pahdxdQ — j  ■ 

dt  (27) 

dNao  dQ2dx 

As  in  the  membrane  case,  we  neglect  the  IVier  order  term  — — .  The  influence  of  static 

<70  2 

tension  will  be  incorporated  as  in  the  e  ?r  membrane  case  by  a  consideration  of  the  deformed 
shape  of  the  shell  and  the  radial  component  that  occurs  due  to  the  deformation  (see  figure  6).  All 
other  dynamic  forces  are  summed  from  an  undeformed  shell  element.  Therefore,  the  T~  term  is 
added  to  the  force  sum  of  equation  (27),  resulting  in 


dx" 


d2w 


„a2w  1  d<2er  dQxr  Nee 

dx2  ade  dx  a  +Pr  P  dt2 


(28) 


Figure  8.  Bending  Shell  Element  FBD  Forces  Only  in  the  r,  x  Plane 

Summing  forces  in  the  longitudinal  x-direction  in  figures  7  and  8  yields 
+ 


2/. 


a2n 

at2 


m  =  pahdxdQ, 


dN 


N xx  +  N xx  +  dx)  °dQ  +  (—  N0x  +  N9x  +  ■^'"'dQ)  dx  +  pxadQdx  =  pahdxdQ 

which  can  be  expressed  as 


BN, 


Qx 


a2« 


ar 


dNxx  ^  1  dN0x 

dx  +  adQ  +P*  =  Ph-- 


d2u 


df 


(29) 
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Figure  9.  Bending  Shed  Element  FBD  Forces  Only  in  the  r,  0  Plane 


Referring  back  to  figures  7  and  9  and  considering  the  forces  acting  in  the  circumferential 
direction,  we  must  now  consider  the  circumferential  component  of  the  shear  force  Q$r: 


XT-  U  V 

=  m57- 


m  =  pahdxdQ, 


3iVnQ  (jQ 

(-^ee  +  ^ee  +  ae  de)dx+  iQ*r  +  dQ 


0N 


X0 


a2v 


+  (-^x9  +  ^jce  +  5_  dx)  ad6  +  peadOdx  =  pahdxdQ—j 
dx  ^ 


which  becomes 


+  1  ^ee  2er 


,a2v 


ajc  T  adO  a 


(30) 
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We  now  consider  the  moments  acting  on  the  shell  element.  Referring  to  figures  7  and  10,  we 
begin  by  summing  moments  about  the  x-  and  0-axes,  setting  the  sum  equal  to  the  mass  moment  of 
inertia  times  the  angular  acceleration  about  the  respective  axis.  The  angular  momentum  is  a  small 
quantity  and  could  have  been  neglected  by  setting  the  moment  sum  equal  to  zero;  however,  to  be 
complete,  it  has  been  included  in  the  derivation.  The  positive  direction  for  a  moment  is  defined  as 
clockwise  about  the  moment  vector  when  we  are  looking  down  the  axis  in  the  positive  direction 
(figure  10).  Beginning  with  the  moments  about  the  jr-axis,  we  have  the  transverse  shear  force  Q$n 
acting  at  a  distance  of  adB/2  from  the  centroid  of  the  shell  element,  along  with  and  Mqq 
yielding  the  following  equation: 

s~\  V  A/f  _  /  mh2 

+  L,Mx  ~  xa^t 2’  lx  =  ’  m  =  P ahdxdQ, 


-  (MxQ  +  ^  dx))  adO  +  (A#ee  -  (Afee  +  ^  ™d6) )  dx 

dQsr  adQ  ph3dxdQd3w 
+  (Ger+(Cer  +  ae  dB))dx—~  =  —2 
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Simplifying  the  above  equation,  we  have 

dx  adQ  +  dr  ~  I2^a^2  (31) 

The  slope  of  the  shell  element  about  the  x-axis  is  expressed  as  and  the  slope  of  the  shell 

O  wv 

element  about  the  0-axis  is 

dx 

About  the  0-axis,  we  have  the  transverse  shear  force  Qxr  acting  through  a  distance  of  dx/2  on 
the  centroid  of  the  element,  along  with  M ^  and  MQx  yielding 

3  2 

=  /e~~2’  7e  =  m  =  pahdxdd, 

A  rHr  i  z 


dQrr  fir 

(-  Qxr  -  (Qxr  +  ^  dx)  )adB™+  (-  M 


0Af 


xx  +  Mxx  +  Yx  dx)adQ 


+  (~mQx  +  mQx 


BMQx 

00 


dB)dx  = 


p  ah3dxdQd3w 

12  0*0t2 


Simplifying  the  above  expression  gives 

_n  dMxxJMQx\ 
Qxr  dx  00  a 


P  A30w 
12  0*0f2 


(32) 


From  equations  (31)  and  (32),  the  transverse  shears  can  be  expressed  in  terms  of  the  moments, 
thus  reducing  the  number  of  stress  resultants  from  ten  to  eight  and  simplifying  the  equations  of 
motion  (equations  (28)  through  (30))  as  follows: 


0/V  xdNdx  d2u 

Tx  +  a00  +Px  =  phtf’ 


^e±l^e0  _l^je  i  0A/ee  ^ _ Jd2v  ,  h 2  d3w 

dx  +a00  aTx  ^00  +Pd  ~  P  i2a2000r2/ 

t,02w_  l^e  l  ^ee  l^e*  nqq _ Jd2w  ,  /j2  04w 

0jc2  a0*0e  a2002  0x2  a002  a  +Pr  P  (dt2  +  12a2deW 


h2  d4w 

12070?/ 


Pr  Pi  Po  Ps  '  (33) 

We  have  not  yet  discussed  pt  or  ps.  Let  us  now  turn  our  attention  to  the  inner  and  outer  fluids 
appearing  in  equations  (24)  and  (33),  which  present  a  loading  on  the  shell.  We  will  derive 
expressions  for  these  pressure  fields  in  terms  of  the  shell  motion  by  applying  boundary  conditions 
at  the  fluid  shell  interface.  Once  we  have  expressions  for  p,  and  ps  in  terms  of  W,  we  will  be  able  to 
solve  for  the  shell  motion  and  subsequently  the  interior  pressure  field.  Lastly,  the  ratio  of  p,{r)IP0 
will  give  us  an  expression  for  the  transfer  function  of  pressure  in  terms  of  wavenumber, 
frequency,  and  the  physical  properties  of  the  shell. 
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FLUIDS  AND  THE  VELOCITY  FIELD  POTENTIAL 

As  the  surface  of  the  shell  is  set  into  motion  by  the  excitation  pressure  wave,  denoted  by  pQ ,  a 
pressure  field,  ps ,  is  generated  in  the  outer  fluid  by  the  shell  surface.  When  the  wavenumber 
associated  with  the  vibration  of  the  shell  surface  is  smaller  than  0)/co,  the  pressure  radiated  into 
the  outer  fluid  is  in  the  form  of  a  propagating  wave.  If  the  wavenumber  of  the  shell  vibration  is 
greater  than  (£>/c0,  the  pressure  field  decays  exponentially  from  the  surface  of  the  shell  in  the  +r 
direction.  A  similar  situation  exists  for  the  inner  fluid  for  p,  based  on  co/c(.  We  will  model  the  fluid 
as  an  ideal  (nonviscous)  linearly  elastic  media  that  cannot  sustain  shear  stress,  even  when  it  is  in 
motion.13  The  equation  of  state  is  then 


\  =  -ptv  ■ 


whose  components  are  given  by  Hooke’s  law  while  we  set  \if=  0: 

ti  i  — 

and 


'11  ~  ‘22  ~  X33  tykk  ~  P 


T12  =  T13  =  X23  =  °‘ 


(34) 


The  following  scalar  velocity  potential  describes  the  dynamic  motion  of  the  fluid14  in  cylindrical 
coordinates: 


where 


V2(p  =  1^/ 

9/  c}*r 

v2  = 

3  r2  r^r  r2dd2  dz2 


(35) 


K  fp 

In  the  case  of  fluids,  J  ^  ,  which  is  equivalent  to  cy  =  J  —  ,  where  (3  is  the  bulk  modulus  of 
the  fluid.  The  particle  velocity  and  the  pressure  in  the  fluid  can  be  expressed  in  terms  of  the 
velocity  potential  as  follows: 


aw 


dt 


f  _ 


=  Vcp, 


(36) 


and 


P  =  ~  P 


a<p, 


fdt 


(37) 
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or 


The  general  solution  to  equation  (35)  is  given  by 
<py  =  Ff(r)  cos  (n0)  e‘^kx~a^ 

where 

Ff(r)  =  DJn{kxr )  +  EYn(klr), 
and 

Ff(r)  =  GIn(k2r)  +HKn(k2r), 


9/  =  ^/(r)  * 


7  CO  7 

=  -2  - k 2 
cf 


-)  0) 
2  _  ,  2  w 
K2  -  K  j 

Cf 


for 


for 


or 


>r, 


*2>^ 

,2 

cf 


(38) 


(39) 


(40) 


At  this  point,  let  us  temporarily  reduce  the  scope  of  our  study  and  consider  only  the  axisymmetric 
zeroth  order  response  ( n  =  0)  in  equations  (38)  through  (40). 


OUTER  FLUID  VELOCITY  FIELD  POTENTIAL 

For  an  outgoing  wave  described  by  e  yk'r  ,  the  field  in  the  outer  fluid  is  physically  a 
radiation  from  the  surface  of  the  cylinder.  The  radial  dependence  F(r)  takes  on  the  form  of  a 
Hankel  function  of  the  first  kind;15  i.e.,  HlJgr)  =  J0(gr)  +  iYJgr).  Therefore,  the  radial  function 
in  the  outer  fluid  is  given  by 

Fs(r)  =  MHl(gr),  g2  =  -  k  for  ^  >*2  • 

C*  C<>  (41) 

2  CO2 

When  k  >  — ,  we  use  the  following  relationship,  substituting  a  modified  Hankel  function  of  real 
argument;16 

Kn(x)  =  ^in+{Hln(ix). 


The  velocity  potential,  <p,  must  approach  0  as  r  — >  oo;  therefore,  10  is  not  an  appropriate  choice 
because  it  approaches  infinity  as  r  — >  <».  Of  the  two  possible  Bessel  functions  in  equation  (40), 
only  the  Ka  Bessel  function  provides  the  necessary  bounded  solution  to  this  problem. 
Consequently,  we  use  the  following  form  for  the  radial  dependence  of  the  velocity  potential: 

.  ,  go2  „  ,  a)2 

Fs(r)  =  HK0(fr),  f  =  k2  -  —  for  *2>—  • 

co  c°  (42) 


The  pressure  component  for  the  radiated  field  is  obtained  with  equations  (37)  and  (38)  while 
we  substitute  an  appropriate  form  for  F(r),  depending  on  the  value  for  k: 

p  (r)  =  /cop  F  (r)ei(kx~w) . 

Fs  y°  sy  '  (43) 

From  equation  (36),  we  can  write  an  expression  for  the  particle  velocity  in  the  fluid: 

8h\ 


-  JL/r  (r\pi{kx~a,) 
dt  "  dr  sK) 


(44) 
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INNER  FLUID  VELOCITY  FIELD  POTENTIAL 

Next,  when  considering  the  behavior  of  the  field  at  r  =  0,  we  make  the  observation  that  Y0  and 
K0  approach  infinity  as  the  argument  of  the  functions  approaches  zero.  However,  because  the 
interior  pressure  field  must  be  finite  at  r  =  0,  these  solutions  are  inappropriate  and  must  be 
discarded.  Bessel  functions  J0  and  IQ  will  be  retained  for  describing  the  velocity  potentials  in  the 
inner  fluid.  Bessel  function  la  replaces  J0  as  the  argument  becomes  imaginary:17 


fn(x)  =  r  Jn  (ix), 


Fj  (r)  =  DJ0(lr), 


Fs(r)  =  Gl 0  (mr) , 


i2  =  V*' 

cf 


0) 


2  a 

m  =  k - - 

cf 


for 


for 


>k 2 

cf 


k2>- 

cf 


(45) 


(46) 


Similar  to  equations  (43)  and  (44),  the  pressure  and  velocity  in  the  inner  fluid  are 

i  ( kx~(Ot ) 


Pi(r)  =  icop.F.(r)e 


^Wi  9  r  /  \  i(kx-(Ot) 

3i  =iF'(r)e 


(47) 


The  derivatives  of  the  Bessel  functions  needed  for  the  fluid  velocity  expressions  of  equations  (44) 
and  (47)  are 


^Asr)  =  ~gH\(gr), 


fK0{fr)  =-fKx(fr), 


±J0{lr)  =  -IJ  x  (Ir) , 

j *-l0(mr)  =  mI\  ( mr )• 


(48) 
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SHELL/FLUID  BOUNDARY  INTERACTION 


AXISYMMETRIC  DEVELOPMENT 


Returning  to  the  shell/fluid  interaction,  we  calculate  the  fluid  pressures  in  terms  of  the  shell 
motion.  This  will  allow  for  the  solution  of  the  shell  displacements,  which  are  now  the  only 
unknowns.  We  will  distinguish  between  inner  and  outer  shell  radii  when  calculating  the  fluid 
velocities.  This  distinction  provides  a  more  accurate  value  for  the  fluid  velocities  and  pressure 
loadings  than  does  the  use  of  the  shell  midradius  for  fluid  calculations.  The  first  step  equates  the 
shell  velocity  with  the  outer  fluid  velocity  at  the  shell  outer  radius,  aQ,  and  thereby  solves  for  the 
constants  M  and  H.  The  second  step  sets  the  shell  velocity  equal  to  the  inner  fluid  velocity  at  r  = 
a,.  Inserting  the  expressions  from  equation  (48)  that  correspond  to  the  outer  fluid  into  equation 
(44)  gives  the  fluid  particle  velocity  as 


dws 

3/ 


Vtps  =  -fHKi(fr)ei(kx-wt) 


(49) 


Differentiating  equation  (26)  to  obtain  the  shell  velocity  in  the  radial  direction  and  then  equating 
the  shell  and  outer  fluid  velocity  at  r  =  a0  while  suppressing  the  exponential  term  yields 

3  wr, 

-i<oW  =  -gMH\  ( ga0 ),  :.M  = 


dw 

3 7 


3  / 


iwW 


gH\  ( ga0 ) 


-/W  =  -fHKx{fa0), 


:.H 


iaW 

/W 


(50) 


In  a  similar  fashion,  we  equate  the  shell  velocity  with  the  inner  fluid  particle  velocity  at  r  =  a,  and 
solve  for  the  constants  D  and  G  as 


3w 

37 


Bws 


-i(oW  =  -IDJX  (laj). 


itoW 


-/(oW  =  mGI^ima^, 


:.G  = 


/0)W 

ml ]  (fflfl;) 


(51) 


The  constants  in  the  expressions  for  the  fluid  pressure  fields  can  be  eliminated,  and  the  following 
simplified  results  are  obtained  using  equations  (41),  (42),  (43),  and  (50)  for  ps  and  equations  (45), 
(46),  (47),  and  (5 1 )  for  p{. 
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Outer  Fluid 


Inner  Fluid 


1  i(kx-m) 


Ps  (r)  = 


-o^p  WW0{gr)e 
gH\  ( ga0 ) 


Mr)  = 


-<a2PiWJ0(lr)el 

Ux  do,) 


i(kx-G)t) 


Ps(r)  = 


-(02poWK0(fr)e‘ 

/* ,  ( fa0) 


i(kx-(Ot) 


Pi(r)  = 


0)2p  WIa  ( mr )  en 
mlx  (mat) 


i(kx-  (0?) 


The  internal  and  external  pressure  fields  have  been  expressed  in  terms  of  the  radial  motion  of  the 
shell,  W.  Let  us  make  the  following  simplifications  to  equation  (52)  and  evaluate  the  pressure  at 
the  required  radii,  then  we  will  return  to  the  shell  equations  of  motion: 


,  .  /  -v  fj7  i(kx-(Dt)  .  .  03  PoHo(8°o)  f_r  ®  v,  fc2 

Ps(ao)  =  «,K) We  as(a0)  =  - — y~ - —  for  ~2  ’ 

gHx{ga0)  cQ 


,  (a  )  .  for  *>> 


, ,  , ,  -®V»('r> 

pi(r)  =  ai(r)H'r  «|(r)  = 


®  ,2 

T  >  A, 

ci 


a  (0  = 


ti>2pJ0(mr) 
ml ]  (ma,) 


k>Za 


,  ,  ,  W./  i(kx-ti)t)  03  P  AoAai)  f 

Pf  ("/)  =  <*AaJWe  «/(«,)  =  -  /yj/-)—  fOT 


^>*2, 


co2P,/0(ri)  fQr 

i{  ,} 
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AXISYMMETRIC  RESPONSE 


MEMBRANE  SHELL 


Let  us  simplify  the  problem  and  consider  only  zeroth  order  axisymmetric  excitation  in  the 
radial  direction.  Pressure  fields  will  be  symmetric  with  respect  to  0;  therefore,  we  will  simplify 
the  equations  with  the  following  conditions  and  eliminate  the  px  excitation: 


=  0,  v  =  0,  pe  =  0,  and  px  =  0  • 

The  specially  orthotropic  stress-strain  relations  are  used,  which  are  represented  by  the  unbarred 
quantities.  Under  these  simplifications,  equations  (22)  through  (24)  reduce  to 


d2u  _  1  dw 
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The  excitation  pressure  field,  pa,  and  the  solution  for  the  displacement  response  are 
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Using  the  derivatives 
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we  write  the  equations  of  motion,  equations  (56)  and  (57),  as 
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Assembling  equations  (60)  and  (61)  in  matrix  form,  we  now  have  the  system  described  in  terms  of 
a  combined  dynamic  stiffness  and  mass  matrix  times  a  displacement  vector  equal  to  a  forcing 
vector: 
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p/jco2  -  Cuk 2 

The  solution  for  W  is 
W  =  - 
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2  ('ll 

By  substitution  of  cn  =  ,  equation  (63)  can  be  rearranged  into  a  form  commonly  seen  in 

some  of  the  literature.  When  the  denominator  is  multiplied  by  1/Cn/1/Cn,  equation  (63)  then 
becomes 
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(64) 


Transfer  Functions  and  Wave  Speeds 

We  can  now  proceed  directly  to  an  expression  for  the  transfer  function  between  the  internal 
pressure  and  the  applied  external  pressure  excitation.  The  pressure  field  in  the  internal  fluid  can  be 
written  directly  from  equations  (54)  and  (64)  as 
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The  magnitude  squared  of  the  transfer  function  can  be  expressed  as 


in2  = 
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P,(r) 
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The  wave  speeds  in  the  system  can  be  analyzed  by  setting  the  determinant  of  the  system 
matrix,  equation  (62),  equal  to  zero  and  rearranging  the  resulting  expression  in  terms  of  c  =  (o/k. 
Even  for  this  relatively  simple  system  of  a  membrane  shell  with  inner  and  outer  fluids,  we  would 
have  to  resort  to  numerical  techniques  to  determine  the  behavior  of  the  wave  speeds. 


We  can  gain  insight  into  the  free  wave  speeds  of  the  structure  by  working  with  the 
denominator  of  equation  (64)  for  a  limited  domain  of  wavenumber  and  frequency.  At  a  resonant 
wavenumber  and  frequency  corresponding  to  the  free  wave  speed  of  the  system,  W  will  reach  a 
maximum  when  the  denominator  equals  zero.  We  will  limit  our  domain  of  interest  to  a 
wavenumber  range  that  corresponds  to  a  small  argument  approximation  for  the  Bessel  functions. 
This  range  is  not  necessarily  “low  wavenumber,”  as  in  the  acoustic  use  of  the  word.  The  region 
will  also  be  one  of  low  frequency.  The  effects  of  the  outer  fluid  and  the  static  tension  will  be 
considered  negligible.  We  will  arrive  at  simple  expressions  for  the  breathing  wave  and 
extensional  wave  speeds. 

2  G)2 

For  the  region  k  >  — ,  we  make  an  approximation  for  p,  by  assuming  that 
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which  results  in 
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We  set  T  =  0  and  consider  p/ico2  to  be  small  compared  with  the  other  terms  because  we 
restrict  ourselves  to  low  frequency.  We  find  the  denominator  of  equation  (64)  simplified  to 
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+  CL.  (a,)  =  0 
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For  wavenumbers  larger  than  klt,  equation  (68)  reduces  to 
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C22  c 


12 


2Pj.c2 


=  0 


a\  Mi '  “/  (69) 

Solving  equation  (69)  for  c,  we  find  the  simplified  expression  for  the  breathing  wave  speed,  cb, 
given  as  a  function  of  the  stiffness  matrix  coefficients  where 


ai 

[c  -CM 

\ 2p  ,.a2 

r22  cj 

cb  = 

*  i  /r\  /*■ —  \  i  «  / 

(70) 

The  simplest  expression  for  the  breathing  wave  speed  is  obtained  with  an  isotropic  shell,  where 
the  material  properties  are  the  same  in  the  three  coordinate  directions.  Such  a  shell  is  termed  an 
unreinforced  shell.  Multiplying  the  elements  of  the  stiffness  matrix,  equation  (5),  by  h  would  then 
give  us  the  stiffness  coefficients  needed  for  equation  (70): 

r  _  hE  hE  „  vhE 

'-'ll  ~ 

(71) 


1  -v 


1  -v 


1  —  v 


For  the  isotropic  shell,  with  all  the  simplifying  assumptions  previously  mentioned,  equation  (70) 
reduces  to  the  following  expression: 


= 


I  hEa-t 

2pp 


(72) 


The  compliance  coefficients  exhibit  a  symmetry  property  that  results  in 


12 


21 

E, 


(73) 


The  specially  orthotropic  stiffness  matrix  in  equation  (4)  simplifies  to  the  following  expression 
when  it  is  substituted  into  equation  (70),  making  use  of  the  symmetry  property  of  equation  (73)  as 
follows 


cb  = 


\hE2ai 

2p,.a2 


(74) 


The  breathing  wave  speeds  given  by  equations  (72)  and  (74)  are  not  affected  by  the  axial 
properties  of  the  shell.  This  condition  is  a  direct  result  of  the  simplifying  assumption  that  k  >  ku. 
In  general,  the  breathing  wave  is  affected  by  the  axial  shell  properties;  however,  the  effect  cannot 
be  observed  from  the  simple  expressions. 


If  the  fluids  and  static  tension  are  ignored  in  the  denominator  of  equation  (64),  the  velocity  for 
the  second  wave  that  the  shell  will  propagate  is  given  by  equation  (75).  Two  terms  with  (o2  have 
been  set  to  zero  since  they  are  small  compared  with  other  terms.  Therefore,  equation  (75)  has 
been  linearized  with  respect  to  frequency.  This  second  wave  is  commonly  referred  to  as  the 
extensional  wave,  because  its  primary  displacement  is  in  the  longitudinal  direction,  representing  a 
stretching  and  contraction  of  the  shell,  as  given  by 
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(75) 
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NONAXISYMMETRIC  RESPONSE 


MEMBRANE  SHELL 
Specially  Orthotropic 

Returning  to  the  full  system  description  seen  in  equations  (22)  through  (26),  we  will  not  make 
the  simplifying  assumptions  that  reduce  the  problem  to  the  zeroth  order  axisymmetric  response. 
In  keeping  the  variation  with  respect  to  0,  we  will  analyze  the  response  of  the  structure  to 
nonaxisymmetric  excitation.  Working  with  the  response  (equation  (26))  allows  the  following 
derivatives  to  be  formed: 

u  =  (/cos  (/i0)  e  v  =  Vsin  (nd)  e  (-kx~(0!\  w  =  W cos  (nQ)  e' ^kx  ~ 
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Equations  (53)  through  (55)  will  be  modified  with  the  addition  of  the  cos(nQ)  term,  as  shown 
in  equation  (38),  along  with  the  generalization  of  the  Bessel  function  to  order  n : 

ps(a0)  =  af(a0)cos  («0)  We  (kx~ (0° 
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The  prime  used  with  the  Bessel  functions  in  equations  (78)  through  (80)  indicates  differentiation 
with  respect  to  r,  according  to  the  following  example 

Equations  (77)  through  (80)  and  equation  (25)  are  substituted  into  equations  (22)  through 
(24).  The  circumferential  shear  stress,  p0,  will  be  assumed  to  be  equal  to  zero.  Retaining  px  allows 
for  the  study  of  the  system  response  to  harmonic  axial  shear  stress  excitation.  We  arrive  at  the 
three  equations  in  the  x,  0,  and  r  directions: 
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Rearranging  equation  (81),  we  have  the  following  system  description  in  matrix  form: 
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Solution  of  equation  (82)  for  W  will  give  us  a  description  of  the  interior  fluid  field  and  the  transfer 
function  between  applied  outer  normal  pressure  and  interior  fluid  pressure  at  radial  position  r, 
which  is 
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The  pressure  in  the  internal  fluid  can  be  written  from  equations  (79)  and  (83)  as 
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Let  us  separate  equation  f84)  into  two  expressions  that  individually  show  the  effect  of  outer 
normal  pressure  on  internal  fluid  pressure  and  the  effect  of  shear  stress  on  internal  fluid  pressure. 
Equation  (85)  is  an  expression  for  the  magnitude  squared  of  the  transfer  function  between 
pressure  in  the  internal  fluid  and  the  applied  normal  pressure  excitation: 
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The  magnitude  squared  of  the  transfer  function  between  pressure  in  the  internal  fluid  and  the 
applied  shear  stress  in  the  longitudinal  jc-direction  is 
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BENDING  SHELL 
Generally  Orthotropic 

Returning  to  the  bending  shell,  we  want  to  solve  equation  (33)  for  the  displacement  magnitude 
W  as  was  done  in  the  last  section  for  the  membrane  shell.  The  excitation  pressure  fields,  pa  and  px, 
and  the  solution  to  equation  (33)  are  given  by  the  following  expressions: 

Po  =  Px  = 


The  interior  pressure  fields  from  equations  (78)  through  (80)  are  modified  accordingly  as 
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The  solution  proceeds  as  follows.  The  displacements,  equation  (87),  are  differentiated  as 
required  and  substituted  into  the  equations  for  the  stress  resultants,  equation  (13).  The  stress 
resultants  are  then  differentiated  as  required  and  substituted  into  the  equations  of  motion,  equation 
(33),  along  with  the  forms  for  the  pressure  fields  p,  and  ps  given  by  equation  (88).  Terms  are 
collected  and  the  following  dynamic  system  matrix  results: 
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The  elements  of  the  system  matrix  S  are 
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We  will  solve  equation  (89)  for  W  twice,  first  using  PQ  as  the  excitation  with  Px  =  0  and  then 
with  Px  as  the  excitation,  with  P0  =  0,  resulting  in 
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The  corresponding  transfer  function,  magnitude  squared,  between  pressure  in  the  inner  fluid 
and  the  applied  external  pressure  P0  or  Px  can  be  written  as  follows: 
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RESULTS 


Let  us  now  apply  the  dynamic  model  derived  thus  far  by  considering  the  behavior  of  specific 
shells.  We  will  examine  the  deformed  shape  of  the  structure  for  the  first  four  modes  of  vibration 
(i.e.,  n  =  0,  1,2,  and  3)  for  the  two  solutions  previously  discussed.  The  first  case  consists  of  a 
standing  wave  in  the  circumferential  0-direction  accompanied  by  a  traveling  wave  in  the 
longitudinal  x-direction;  the  second  case  consists  of  a  traveling  wave  solution  in  0  and  x. 

Next,  using  the  solutions  represented  by  equations  (92)  and  (93),  we  shall  examine  the  system 
response  as  a  function  of  longitudinal  wavenumber  and  frequency  for  each  elastic  model  at  each 
of  the  first  four  modes  of  propagation.  Effects  of  fluid  loading  and  reinforcement  wrap  angle  <|> 
will  be  explored  as  well  as  the  effects  of  an  internal  core  for  the  n  =  0  mode. 

We  will  then  explore  the  strain  in  an  air-backed  shell  as  developed  by  the  normal  pressure 
excitation.  With  this  structure,  we  shall  increase  the  stiffness  of  the  shell  and  observe  the  effects 
on  the  magnitude  of  the  strain  response  as  a  function  of  longitudinal  wavenumber  and  frequency. 

MODE  SHAPES 

The  mode  number  n  is  defined  as  the  number  of  wavelengths  of  oscillation  that  a  quantity 
undergoes  with  respect  to  the  circumferential  direction  at  a  constant  x-location.  Figure  1 1  depicts 
how  the  displacement  u,  v,  or  w  would  vary  with  respect  to  the  circumference  of  the  shell.  We  note 
that  the  n  =  0  shape  is  such  that  all  points  on  the  circumference  move  in  phase  with  one  another 
for  a  given  x-location.  For  n  =  1,  we  observe  one  wavelength  around  the  circumference.  As  the 
mode  number  increases  (i.e.,  n  =  1,  2,  3...),  more  wavelengths  are  added  within  the 
circumferential  length  n d.  In  this  report,  two  solutions  to  the  equations  of  motion  for  the  shell  and 


fluid  are  given  by  equations  (26)  and  (87).  Both  equations  have  been  applied  to  depict  the 
deformation  of  a  shell  at  a  given  instance  in  time.  Remember  that  these  solutions  are  for  an 
infinite  cylinder;  therefore  the  mode  shape  are  constructed  with  a  length  of  shell  equal  to  two 
wavelengths  in  order  to  visualize  the  displacement  that  the  shell  is  undergoing.  This  length  is 
arbitrarily  chosen  at  two  wavelengths.  Although  another  length  would  have  been  equally  correct, 
it  would  not  necessarily  be  as  instructive  for  visualizing  the  mode  shapes. 

Standing  Wave 

Figures  12  through  28  have  been  created  from  equation  (26)  for  the  displacements.  As  stated 
before,  equation  (26)  represents  a  standing  wave  solution  in  0.  The  mode  shape  results  from  a 
wave  traveling  in  the  clockwise  and  counterclockwise  directions.  Referring  to  figure  13,  the  net 
result  of  these  pressure  fields  on  the  shell  is  a  vibration  occurring  in  the  x,  0  =  0  plane.  The  w 
component  of  mode  1  causes  a  translation  of  the  circumference  of  the  shell  from  the  center, 
resulting  in  a  flexural  mode  of  vibration. 

Figure  12  depicts  the  classic  breathing  wave,  as  all  points  around  the  circumference  at  a  given 
jc  location  move  in  phase  in  the  radial  direction.  Figures  14  and  15  depict  two  and  then  three 
wavelengths  of  variation  in  w  around  the  circumference.  The  extensional  wave,  which  is  the  u 
component  of  the  n  -  0  mode  of  propagation,  is  shown  in  figure  16,  with  the  higher  order  modes 
of  u  shown  in  figures  17  through  19.  The  zeroth  order  mode  in  v  results  in  zero  tangential 
displacement  as  seen  from  equation  (26)  and  figure  20.  Figures  21  through  23  show  the  effects  of 
the  higher  modes  of  v  on  the  shell. 

The  effects  of  all  three  components  of  displacement  (w,  u,  and  v)  are  imposed  on  the  shell  for 
each  mode  in  figures  24  through  27.  To  visualize  the  effects  of  the  combined  displacements,  the 
magnitudes  of  w,  u,  and  v  have  been  arbitrarily  chosen  to  be  equal  to  one  another.  Actually,  the 
amplitudes  U  and  V  are  much  smaller  than  W  for  a  radial  pressure  excitation,  which  we  will 
consider  later.  If  the  actual  amplitudes  had  been  displayed  for  the  radial  pressure  excitation  in  the 
correct  proportion  in  figures  24  through  27,  only  the  w  mode  shape  would  be  visible.  Furthermore, 
the  exact  proportion  of  the  displacements  W/W,  V/W,  and  U/W  are  functions  of  the  shell  and  fluid 
and  will  be  different  for  each  shell  model  (i.e.,  isotropic,  specially  orthotropic,  generally 
orthotropic).  These  proportions  also  depend  on  the  magnitude  of  the  elastic  constants.  We 
emphasize,  therefore,  that  figures  24  through  27  provide  qualitative,  not  quantitative,  information. 

Figure  28  is  an  animation  of  the  w  component  of  the  n  =  2  mode,  showing  the  mode  shape  for 
six  snapshots  in  time  for  the  wave  propagating  in  the  +;c-direction.  Careful  inspection  of  the  figure 
indicates  that  the  displacement  wave  is  propagating  to  the  left,  or  positive  x-direction. 
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Mode  1 

Radial  Component  w 


Figure  13. 6  Standing  Wave  Mode  Shape;  n  =  1,  w  Component 
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Figure  14. 6  Standing  Wave  Mode  Shape;  n  =  2,w  Component 


Figure  15. 6  Standing  Wave  Mode  Shape;  n  =  3,w  Component 
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Mode  0 

Longitudinal  Component  u 


T 


Figure  16. 0  Standing  Wave  Mode  Shape;  n  =  0,  u  Component 


Figure  17. 6  Standing  Wave  Mode  Shape;  n  =  l,u  Component 
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Figure  18. 6  Standing  Wave  Mode  Shape;  n  =  2,u  Component 


Mode  3 

Longitudinal  Component  u 


Figure  19. 6  Standing  Wave  Mode  Shape;  n  =  3,u  Component 
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Mode  0 

Tangential  Component  v 


Mode  1 


Figure  21. 6  Standing  Wave  Mode  Shape;  n  =  1,  v  Component 
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Mode  3 

Tangential  Component  v 


Figure  23. 6  Standing  Wave  Mode  Shape;  n  =  3,  v  Component 
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Figure  24. 0  Standing  Wave  Mode  Shape;  n  =  0  for  w,  u,  and  v  Components 


Figure  25. 0  Standing  Wave  Mode  Shape;  n-  1  for  w,  u,  and  v  Components 
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Figure  26. 6  Standing  Wave  Mode  Shape;  it  =  2  for  w,  u,  and  v  Components 


Mode  3 

Components  w,u,&v 


Figure  27. 6  Standing  Wave  Mode  Shape;  n  =  3  for  w,  u,  and  v  Components 
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Figure  28.  Animation  of  6  Standing  Wave  Mode  Shape;  n  =  2,w  Component 


Traveling  Wave 

Similar  to  the  previous  section,  figures  29  through  45  have  been  created,  from  equation  (87) 
and  depict  mode  shapes  for  the  first  four  modes  of  a  traveling  wave  excitation  in  x  and  0.  The 
direction  of  propagation  is  in  the  +x-  and  +0-directions. 

It  is  noteworthy  that  the  mode  shapes  spiral  around  the  circumference  of  the  shell  in  the  +0- 
direction  as  the  wave  progresses  in  the  +jc-direction.  The  w  component  of  mode  1  becomes  a 
bending  vibration  that  causes  a  corkscrewing  of  the  shell  rather  than  an  in-plane  vibration  as 
occurs  in  figure  13  for  the  standing  wave  mode  shape. 

Special  attention  is  warranted  for  the  v  component  of  mode  n  =  0  seen  in  figure  37.  As  was 
stated  in  the  section  on  Axisymmetric  Response  Membrane  Shell,  v  =  0;  this  condition  is  true  for 
isotropic  and  specially  orthotropic  formulations  for  the  shell.  The  only  time  that  v  is  not  equal  to 
zero  for  the  zeroth  mode  n=0  is  in  the  generally  orthotropic  case.  Therefore,  figure  37 
corresponds  to  a  generally  orthotropic  shell  only.  Isotropic  and  specially  orthotropic  shells  are 
depicted  by  figure  20  for  the  traveling  wave  case. 

The  combined  effects  of  the  three  components  of  the  displacement  are  shown  in  figures  41 
through  44  (where  the  displacement  magnitudes  are  equal,  as  in  the  standing  wave  case)  only  to 
allow  for  visualization  of  the  effects  of  displacement.  The  traveling  wave  animation  sequence  is 
displayed  in  figure  45;  the  wave  is  traveling  in  the  +x-  and  +0-direcfi<  »ns. 
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Figure  29. 0  Traveling  Wave  Mode  Shape;  n  =  0,  w  Component 


Figure  30. 0  Traveling  Wave  Mode  Shape;  n  =  l,w  Component 
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Mode  0 

Longitudinal  Component  u 


Figure  33. 6  Traveling  Wave  Mode  Shape;  n  =  0,«  Component 


Mode  1 

Longitudinal  Comnonent  u 


Figure  34. 6  Traveling  Wave  Mode  Shape;  n  =  1,  u  Component 


55 


Mode  2 

Longitudinal  Component  u 


Mode  3 

Longitudinal  Component  u 


Figure  36. 0  TVaveling  Wave  Mode  Shape;  n-\u  Component 


56 


Mode  0 

Tangential  Component  v 


Figure  38. 6  Traveling  Wave  Mode  Shape;  n-  1,  v  Component 
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Figure  39. 0  Traveling  Wave  Mode  Shape;  n  =  2,  v  Component 


Mode  3 

Tangential  Component  v 


Figure  40. 6  Traveling  Wave  Mode  Shape;  n  =  3,  v  Component 
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Figure  41.  6  Traveling  Wave  Mode  Shape;  n  =  0  for  w,  u,  and  v  Components 


Mode  1 


Figure  42.  6  Traveling  Wave  Mode  Shape;  n  -  l  for  w,  u,  and  v  Components 


Figure  43. 6  Traveling  Wave  Mode  Shape;  n  =  2  for  w,  u,  and  v  Components 


Mode  3 

Component  w,u,&v 


Figure  44. 6  Traveling  Wave  Mode  Shape;  n  =  3  for  w,  u,  and  v  Components 
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Figure  45.  Animation  of  6  Traveling  Wave  Mode  Shape;  n  =  2,w  Component 


MAGNITUDE  OF  THE  SHELL  RESPONSE 


Now  that  we  have  examined  the  deformation  of  the  shell,  let  us  consider  the  magnitude  ( W ’  U, 
and  V)  of  the  shell  response,  normalized  by  the  excitation  pressure  or  stress,  as  a  function  of 
wavenumber  and  frequency. 

Urethane  Shell  Parameters 

We  will  proceed  by  using  material  properties  that  approximate  an  ordinary  garden  hose  for  the 
first  set  of  numerical  results.  In  subsequent  sets,  only  the  parameters  that  are  changed  will  be 
described  as  the  other  parameters  will  remain  the  same. 

Matrix 

E  =  1.5  x  107  Pa  C  =0.01 

v  =  0.45  dimensionless 

m 

pv  =  1000  kg/m3 

Reinforcing  Fiber 

Ef=  1.5  xlO9  Pa  ^=0.0 

=  0.15  dimensionless 

0.9907  x  10-3  meters  (0.039  inches)  outer  diameter  each  strand 
12  strands 

Shell 

a  =  0.01270  meters  (0.5  inches) 
h  =  0.00317  meters  (0.125  inches) 

Fluids 

p,  =  1 000  kg/m3  Cj  =  1 500  m/s 

p0  =  1000  kg/m3  ca  =  1500  m/s 
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Structural  Damping 


We  have  not  yet  discussed  loss  within  the  structure.  The  well-known  structural  damping 
model  will  best  describe  the  dissipation  of  energy  within  the  shell.  Structural  damping  introduces 
a  complex  modulus  of  elasticity  where  Em  and  Ef  become 

E  —  E  ( 1  +  C  / ) , 

m  m v  ' 


E/-£/<1+V')- 

The  loss  tangent  or  structural  loss  factor,  £,  is  defined  as  the  ratio  of  the  in-phase  to  out-of-phase 
components  of  the  modulus.  The  loss  factor  is  determined  experimentally  for  the  material  of 
interest.  Such  experiments  usually  involve  excitation  of  a  specimen  with  a  harmonically  applied 
force  and  measurement  of  the  in-phase  and  out-of-phase  displacement  responses  with  respect  to 
the  applied  force. 

Based  on  the  previously  described  material  properties,  the  rule  of  mixtures  has  been  applied 
(appendix  B)  to  calculate  the  composite  properties  for  the  reinforced  shell.  Where  a  plot  or  image 
is  labelled  isotropic,  the  calculation  is  performed  without  reinforcement.  When  labeled  specially 
orthotropic,  the  reinforcing  fibers  are  included.  The  angle  <j)  will  always  be  0  for  the  isotropic  and 
specially  orthotropic  case.  The  generally  orthotropic  case  has  an  angle  that  is  formed  between  the 
reinforcement  and  the  longitudinal  axis  of  the  shell  (<|>),  as  seen  in  figure  3. 

Displacement  Response 

In  figure  46,  displacements  W,  U,  and  V  have  been  obtained  from  equation  (89),  which 
illustrates  an  isotropic  case,  normalized  by  Pa  with  Px  =  0  at  k  =  58.3  rad/m.  At  the  first  resonance 
(210  Hz),  the  amplitude  of  W  is  the  largest  among  the  three  amplitude  components  W,  U,  and  V. 
There  are  corresponding  displacements  in  U  and  V”  due  to  the  210-Hz  resonance  in  W,  but  they  are 
lower  in  amplitude.  The  resonance,  therefore,  corresponds  to  the  w  component  of  the  n  =  2  mode 
of  propagation.  The  W  response  couples  into  the  U  and  V  displacements,  but  at  a  lower  amplitude. 
Similar  behavior  occurs  at  the  second  resonance,  2000  Hz,  which  corresponds  to  the  u  component 
of  the  n  =  2  mode  of  propagation.  Coupled  responses  are  found  in  the  W  and  V  displacements  as 
well.  The  v  component  of  mode  n  =  2  occurs  at  approximately  4000  Hz,  causing  both  U  and  W 
coupled  displacement  responses. 

The  antiresonances  are  caused  by  a  difference  in  phase  between  the  displacement  of  a  given 
component  of  the  n  =  2  mode  (w  for  example)  and  the  coupling  from  the  u  and  v  components  of 
the  n  =  2  mode  into  w.  In  figure  46,  two  antiresonances  observed  in  the  W  response  occur  under 
the  U  and  V  resonant  peaks,  at  2000  and  4000  Hz,  respectively. 

The  frequencies  and  wavenumbers  corresponding  to  the  peak  values  of  each  component  of  the 
n  =  2  mode  are  members  of  branches.  Branches  will  be  discussed  in  a  later  section  where  we 
observe  a  surface  formed  from  the  W  response.  We  will  be  able  to  discuss  all  three  branches  of  the 
n  =  2  mode  from  considering  a  single  transfer  surface  due  to  coupling. 
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Now  let  us  turn  our  attention  to  the  W  response  in  the  100-  to  1500-Hz  region  (figure  46).  This 
region  can  be  divided  into  three  parts,  from  low  to  high  frequency.  The  first  region  is  a  flat, 
straight-line  stiffness  controlled  region,  which  would  become  increasingly  flat  if  the  next  decade 
below  100  Hz  was  displayed.  The  second  region  is  marked  by  the  resonance  with  increasing 
amplitude  and  phase  difference.  The  third  is  a  straight,  but  sloping,  region  where  the  response  of 
the  system  is  controlled  by  its  mass. 


Figure  46.  Three  Components  of  Mode  n  =  2  at  k  =  58.3  rad/m 


An  interesting  comparison  to  the  field  of  modal  analysis  can  be  made  here.  Individually,  the 
modes  of  a  receptance  frequency  response  function  (FRF),  displayed  as  a  log,  log  plot,  are  subject 
to  the  same  three  regions  applied  to  figure  46.  Furthermore,  the  antiresonances  appear  in  the 
receptance  FRF  due  to  a  phase  difference  between  adjacent  modes.  Mode  numbering  is  different 
in  modal  analysis.  In  modal  analysis,  the  resonances,  or  modes,  are  numbered  with  respect  to 
increasing  frequency,  apart  from  a  consideration  of  the  mode  of  wave  propagation.  Therefore, 
several  modes  of  wave  propagation  could  be  supported  by  a  structure.  The  subsequent  numbering 
of  modes  in  the  modal  analysis  would  not  necessarily  correspond  to  the  numbering  of  the  modes 
of  wave  propagation  as  discussed  in  this  work.  Figure  46  depicts  the  situation  for  a  mode  of  wave 
propagation,  (n  =  2)  that  the  structure  and  fluids  will  support. 

From  figure  46,  the  frequencies  of  vibration  for  the  various  components  of  the  mode  shapes 
can  be  determined  from  the  resonant  peaks.  The  correct  displacement  amplitudes  are  also 
available  in  this  figure.  We  see  from  the  first  resonance  in  W  that  corresponding  displacement 
magnitudes  in  U  and  V  are  some  3  and  10  dB  lower,  respectively.  Therefore,  at  approximately  220 
Hz,  which  corresponds  to  the  w  component  of  mode  n  =  2,  we  would  also  observe  the  u  and  v 
components  of  mode  n  =  2  occurring  with  appropriately  reduced  amplitudes  of  vibration. 
Simultaneously  seen  on  the  structure,  but  at  2000  Hz,  is  the  U  resonance,  but  it  is  some  55  dB 
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lower  than  the  value  of  the  amplitude  for  the  resonance  in  W. 

We  see  that  the  amplitude  ratios  depend  on  the  physical  properties  of  the  structure.  This  is 
why  the  displacement  amplitudes  were  arbitrarily  chosen  to  be  equal  in  displaying  the  mode 
shapes;  several  of  the  modes  would  not  be  observable  if  the  amplitudes  were  unequal. 

TRANSFER  SURFACES  IN  THE  k, co  PLANE 

When  equation  (92)  is  evaluated  over  a  range  of  wavenumbers  and  frequencies,  we  have  the 
wavenumber-frequency  plane  shown  in  figure  47,  where  the  color  white  indicates  maximum 
amplitude  and  black  indicates  minimum  amplitude. 
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Figure  47,  Bending  Shell  Transfer  Surface  (Image);  dB  =  101og((/’((«/2)//’o)2) 
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The  image  shown  in  figure  47  is  illustrated  as  a  wire  frame  surface  in  figure  48  to  help  the 
reader  visualize  the  three-dimensional  nature  of  the  response.  The  transfer  surfaces  in  the 
remaining  figures  will  be  shown  as  the  grey  scale  image  seen  in  figure  47.  The  reader  should  refer 
to  this  pair  of  figures  to  interpret  surfaces  in  later  figures. 

Branches 

In  figure  47,  we  see  three  approximately  v,  or  horseshoe-shaped,  curves  or  ridges  symmetric 
about  k  =  0.  The  peak  values  along  a  given  ridge  line  make  up  the  collection  of  points  that 
comprise  a  branch.  Both  the  positive  and  negative  wavenumber  halves  will  be  considered 
members  of  the  same  branch.  Therefore,  any  reference  made  to  a  branch  means  both  halves.  A 
branch  displays  the  relationship  between  the  wavenumber  and  frequency  for  a  particular  mode  of 
propagation.  Thus,  from  left  to  right,  the  branches  are  associated  with  the  w,  «,  and  v  components 
of  the  mode  of  propagation  n  -  1.  Mode  shapes  shown  in  figures  13,  17,  and  21  correspond  to 
these  three  branches,  respectively.  Strictly  speaking,  a  branch  is  one  set  of  roots  belonging  to  the 
determinant  of  the  dynamic  system  matrix  of  equation  (89).  The  tranrfer  surfaces  have  been  taken 
one  step  further  than  just  a  display  of  the  roots  of  the  equation.  Here,  the  magnitude  of  the 
displacement  has  been  calculated,  as  mentioned  earlier.  Velocity  information  is  still,  nonetheless, 
obtainable  from  the  transfer  surfaces  when  we  consider  the  peak  values  of  the  transfer  surface 
ridges  as  representing  the  collection  of  points  that  comprise  the  branches.  When  reference  is  made 
to  branches,  henceforth,  the  peak  value  of  the  ridges  is  meant. 

Phase  Velocity 

The  inverse  of  the  slope  calculated  by  use  of  the  two  points  consisting  of  the  origin  of  the  k, to 
plane  and  any  point  on  a  branch  is  called  the  phase  velocity.  If  the  phase  velocity  does  not  depend 
on  wavenumber,  but  is  rather  a  constant,  the  branches  are  given  by  straight  lines  in  the  k,(&  plane 
and  are  termed  nondispersive.  If  the  phase  velocity  is  a  function  of  wavenumber,  the  branch 
describes  a  dispersive  form  of  wave  propagation. 

Group  Velocity 

The  quantity  ^  calculated  at  a  point  on  a  branch  is  called  the  group  velocity  for  the  particular 

U  ri 

point.  The  group  velocity  defines  the  rate  at  which  energy  propagates  in  the  particular  mode.  For 
nondispersive  modes,  the  phase  and  group  velocities  are  equal.  Therefore,  energy  propagates  at  a 
constant  velocity  for  all  wavenumber  and  frequency  pairs  along  the  nondispersive  branch. 

All  three  branches  in  figure  47  are  dispersive.  The  least  dispersion  is  exhibited  by  the  first 
branch,  which  is  the  farthest  left.  This  branch  actually  becomes  nondispersive  for  sufficiently  high 
wavenumbers. 

Cutoff  Frequency 

Cylindrical  shells  have  a  finite  circumferential  length.  It  is  the  restriction  of  a  spatial 
dimension  that  causes  the  shell  to  become  a  waveguide.  Waveguides  exhibit  cutoff  frequencies 
below  which  particular  modes  will  not  propagate.  This  behavior  is  most  clearly  seen  in  figure  47 
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for  the  second  and  third  branches  of  the  n  =  1  mode.  A  layer  of  finite  thickness  is  also  a 
waveguide  that  exhibits  similar  behavior. 

At  the  cutoff  frequency,  the  group  velocity  equals  zero,  and  energy  does  not  propagate  along 
the  longitudinal  axis  of  the  shell.  Along  the  branch,  as  the  frequency  is  increased  beyond  the 
cutoff  frequency,  the  group  velocity  increases  from  zero,  and  energy  begins  to  propagate  down  the 
longitudinal  axis  of  the  shell. 

We  have  been  considering  propagating  excitations  and  responses  of  the  form 


cos  (/?0)  e 


i(kx-(ot) 


e 


i  ( k.x  +  ;i0  -  (i)l) 


(94) 


When  a  branch  displays  a  cutoff  frequency,  we  find  that  the  behavior  of  the  structure  can  be 
observed  below  cutoff  frequency  by  considering  imaginary  wavenumbers.  Inserting  ik  into 
equation  (94)  gives 


cos (nd) e 


-k.x ( «0  -  to/) 


(95) 


The  imaginary  wavenumber  turns  equation  (95)  into  an  inhomogeneous  vibration  whose 
amplitude  is  decaying  with  increasing  a.  In  this  case,  k  is  actually  an  attenuation  factor,  not  a 
spatial  frequency  as  in  the  propagating  wave  case.  In  the  0-direction  of  equation  (95),  we  first 
have  a  standing  wave  and  then  a  traveling  wave.  The  amplitude  of  the  response  is  shown  in 
figures  49  and  50  for  a  shell  only.  The  effects  of  the  fluids  have  been  removed. 

In  further  discussions,  we  will  explore  the  effects  of  various  parameters  on  the  cutoff 
frequencies  as  well  as  on  the  general  behavior  of  the  transfer  surfaces. 

Isotropic  n  =  2  0  0 


2000  1000 
Frequency  Hz 

Figure  49.  Transfer  Surface  (Real  Wavenumber);  dB  =  10log((W/f>0)2) 
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2000  4000 


Frequency  Hz 

Figure  50.  Transfer  Surface  (Imaginary  Wavenumber);  dB  =  101og((W//*o)2) 

RADIAL  PRESSURE  TRANSFER  SURFACES 
Membrane  Shell 

Figures  51  through  55  describe  the  membrane  shell  dynamic  response.  Figure  52(a)  contains 
two  branches.  The  major  and  lowest  speed  branch  corresponds  to  the  w  component  of  mode  n  =  0. 
whose  shape  is  displayed  in  figure  12.  This  branch  describes  a  mode  of  propagation  referred  to  as 
the  breathing  wave.  The  u  component  of  mode  n  =  0.  whose  displacement  shape  is  seen  in  figure 
16,  is  faintly  visible  in  figure  52(a)  at  a  higher  phase  velocity.  This  u  component  is  referred  to  as 
an  extensional  wave.  From  figure  52(a),  it  is  clear  that  the  radial  P„  excitation  does  not  generate  a 
significant  extensional  wave  level. 

Observing  mode  n  =  1,  we  can  see  the  effect  of  longitudinal  stiffness  on  all  three  branches, 
most  visibly  on  the  second  and  third.  The  group  velocity  is  increased  for  the  third  branch  and  over 
most  of  the  second  branch  as  well  (figure  53(b)). 

The  branches  of  the  generally  orthotropic  surfaces  exhibit  very  interesting  behavior.  For  the 
first  time,  we  observe  an  asymmetry  about  k  =  0.  Figure  51  depicts  the  resultant  propagation 
directions  that  correspond  to  the  plus  and  minus  wavenumber  half  planes.  For  the  plus 
wavenumber  half  plane,  the  +.v,  +9  propagation  direction  rotates  into  the  direction  aligned  with 
the  £2  material  property.  Arriving  from  the  -  v.  +0  direction  (the  negative  wavenumber  half  plane), 
the  propagation  is  principally  in  the  E\  material  property  direction.  For  the  generally  orthotropic 
case,  £-,  *  £  ( ;  therefore,  different  stiffness  is  presented  by  the  structure  to  the  excitation,  and  an 
asymmetry  results.  It  is  interesting  to  note  that  changing  the  sign  of  n  from  positive  to  negative 
will  flip  the  surface  180  degrees  about  the  k  =  0  axis,  as  would  be  expected  from  the  preceding 
discussion. 

Looking  at  mode  n  =  2,  figure  54(a)  -(c).  we  see  that  the  first  branch  does  not  exhibit  a  cutoff 
frequency.  This  behavior  will  be  contrasted  in  the  next  section  on  the  bending  shell,  where  the 
first  branch  of  the  second  mode  exhibits  a  cutoff  frequency. 
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In  general,  for  a  given  branch  of  the  membrane  shell  and  shell  stiffness  type  (i.e.,  isotropic), 
the  cutoff  frequency  increases  with  increasing  mode  number,  and  for  a  given  wavenumber,  the 
group  velocity  slows  down  with  increasing  mode  number. 


Figure  51.  Generally  Orthotropic  Shell  Schematic,  Propagation  Directions 


Bending  Shell 

Figures  56(a)  through  62  correspond  to  the  bending  shell  dynamic  response.  Two  major 
differences  are  observed  in  the  bending  shell  surfaces.  The  first  concerns  the  breathing  wave, 
which  is  the  w  component  of  mode  n  =  0  (see  figure  56(a)-(c)).  At  high  wavenumbers,  the  group 
velocity  increases,  in  contrast  to  the  membrane  case  where  group  velocity  slows  with  increasing 
wavenumber. 

The  w  component  of  mode  n  =  2.  exhibits  a  cutoff  frequency  (see  figure  58(a)-(c)).  We 
observe,  therefore,  the  importance  of  bending  stiffness  on  this  particular  branch. 

As  in  the  membrane  case,  cutoff  frequency  increases  with  increasing  mode  number  for  a  given 
branch.  For  the  u  and  v  components  of  the  modes,  we  see  the  group  velocity  decreasing  with 
increasing  mode  number  for  all  dispersive  branches. 
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Membrane  Shell  Sequence  Continued 

Isotropic 


Frequency  Hz 
(a) 

S.  Orthotropic  n  = 


Frequency  Hz 

(b) 

G.  Orthotropic  n  = 

=  2 

Frequency  Hz 
(C) 

Figure  53.  Membrane  Shell  Transfer  Surface  (n  =  1);  dB  =  10Iog((/,,<a/2)//,o)2) 


Frequency  Hz 
(c) 

Figure  54.  Membrane  Shell  Transfer  Surface  (n  =  2);  dB  =  10log((PI<a/2)//,o)2) 


Membrane  Shell  Sequence  Continued 
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Figure  56.  Bending  Shell  Transfer  Surface  (n  =  0);  dB  =  101og((P,<a/2)/Fo)2) 
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Bending  Shell  Sequence  Continued 
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Figure  57.  Bending  Shell  Transfer  Surface  (n  =  1);  dB  =  10log((P,<a/2)/Fo)2) 
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Bending  Shell  Sequence  Continued 


G.  Orthotropic  n=  3  0=  20 
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Figure  59.  Bending  Shell  Transfer  Surface  (n  =  3);  dB  =  10log((P,{a/2)/Po)2) 


Outer  Fluid  Effects  on  the  w  Component  of  the  n  =  2  Mode 


Figure  62.  Comparison  of  Levels  at  k  =  0  From  Figures  60  and  61 


Effects  of  Wrap  Angle  0  on  the  n  =  2  Mode 

Figure  63  displays  the  effects  of  increasing  the  wrap  angle  from  0  to  90  degrees.  The 
asymmetry  introduced  by  this  angle  disappears  again  at  <j)  =  90  degrees,  where  the  structure 
becomes  once  again  specially  orthotropic. 

The  cutoff  frequency  for  the  first  branch  has  increased  at  <j)  =  90  degrees  because  of  the 
increased  stiffness  in  the  circumferential  direction.  The  third  branch,  the  v  component  of  mode  n 
-  2,  is  moving  in  the  direction  of  increasing  frequency  for  the  same  reason  and  actually  moves  out 
of  the  range  displayed  in  the  k, to  plane. 
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rad/m  k  rad/m 


Effects  of  Wrap  Angle  <j)  on  the  n-2  Mode  Continued 


G.  Orthotropic  n=  2 
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G.  Orthotropic  n= 
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S.  Orthotropic  n  = 


Frequency  Hz 


Figure  63.  Bending  Shell  Transfer  Surface  (n  =  2)  for  0  <  4>  <  90;  dB  =  lOlog ((P,{a/2)/P0)2) 

(Cont'd) 
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Interior  Partial  Cross  Section  Filling  Core 

In  figures  64  and  65,  we  see  the  branch  corresponding  to  the  w  component  of  mode  n  =  0. 
Since  the  branches  are  dispersive,  we  will  describe  their  velocity  by  referring  to  the  phase  velocity 
at  45  rad/m.  The  phase  velocity  in  figure  64  is  37  m/s  and,  in  this  case,  the  shell  is  fluid  filled 
without  an  internal  core.  Let  us  now  place  a  compressible  infinitely  long  core  into  the  shell;  refer 
to  figure  C-l  in  appendix  C.  The  core  has  the  following  parameters: 

r(,  =  7.62  x  lCf  3  m  (0.3  inches) 

v  =  0.45, 

E  =  1.5  x  107  Pa,  C  =  0.01. 

The  transfer  surface  of  figure  65  is  computed  by  substitution  of  equations  C-l  1  and  C-l 2  for 
equations  (54)  and  (55)  when  we  evaluate  equation  (92)  at  n  =  0.  As  seen  in  the  development  in 
appendix  C,  and  now  from  figure  65,  the  effect  of  the  core  is  to  slow  the  phase  velocity.  The  phase 
velocity  from  figure  65  at  45  rad/m  is  28  m/s. 


Frequency  Hz 

Figure  64.  Bending  Shell  Transfer  Surface  Without  Core;  dB  =  10log((P,-(Jo/4)//>o)2) 
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rad/ m 


Figure  65.  Bending  Shell  Transfer  Surface  With  Core;  dB  =  10log((P,<3a/4)/Po)2) 
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An  extreme  case  of  a  compliant  core  is  represented  in  figure  66.  Material  properties  and 
dimensions  of  this  core  are 

rc  =  7.62  x  10-3  m  (0.3  inches), 

v  =  0.45, 

C 

E  =  1.0  x  105  Pa,  C  =  0.01. 

This  core  is  similar  to  soft  foam  rubber,  where  the  fluid  does  not  penetrate  the  foam.  In  this 
case,  the  phase  velocity  has  slowed  to  16  m/s.  This  very  soft  core  gives  rise  to  branches  that  bear 
a  close  resemblance  to  the  air-backed  case  that  we  shall  study  in  a  later  section. 


200  400  600  800  1000 


Frequency  Hz 


Figure  66.  Bending  Shell  Transfer  Surface  With  Soft  Core;  dB  =  10log((P,-(5a/4)/Po)2) 
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AXIAL  SHEAR  STRESS  TRANSFER  SURFACES 


Bending  Shell 

In  this  section,  we  will  explore  the  behavior  of  the  branches  for  the  first  four  modes  of 
propagation  from  the  transfer  surfaces  computed  using  equation  (93).  In  this  case,  the  shell  is 
excited  with  a  longitudinal  shear  stress. 

The  branches  observed  in  figure  67(a)  correspond,  as  before,  to  the  vr  and  u  components  of  the 
n  =  0  mode  of  propagation;  the  first  (slower  speed)  branch  represents  the  breathing  wave  and  the 
second  (faster)  represents  the  extensional  wave.  Figure  67(a)-(c)  should  be  compared  with  figure 
56(a)-(c)  to  note  the  difference  in  structural  response  between  an  axial  shear  stress  excitation  and 
a  normal  pressure  excitation.  In  figure  56(a)-(c),  the  extensional  wave  is  barely  excited,  existing  at 
a  very  low  level.  However,  in  figure  67(a)-(c),  the  extensional  wave  is  strongly  evident.  The 
breathing  wave  is  evident  as  well,  but  it  is  lower  in  level. 

Of  all  the  modes  of  propagation  considered  for  the  shell,  the  extensional  wave  is  the  most 
nondispersive,  mainly  because  this  mode  engages  the  shell  in  minimum  bending  and  very  little 
interaction  with  the  fluids.  The  addition  of  longitudinal  reinforcement  has  increased  the  phase 
velocity  from  138  m/s  in  the  isotropic  case  of  figure  67(a)  to  295  m/s  in  the  specially  orthotropic 
case  of  figure  67(b).  Phase  and  group  velocities  of  the  breathing  wave  have  increased  most 
noticeably  for  the  high  wavenumbers. 

It  is  interesting  to  note  the  generally  orthotropic  case  of  figure  67(c),  where  another  branch  has 
appeared.  The  higher  modes  (n  =  1  through  n  -  3)  shown  in  figures  68  through  70,  display  the 
mode  number,  cutoff  frequency,  and  group  velocity  characteristics  that  their  counterparts  did  in 
figures  57  through  59. 

Effects  of  Wrap  Angle  <J>  on  the  n-0  Mode 

Let  us  consider  mode  n  =  0  and  explore  the  effects  of  varying  the  angle  of  the  reinforcement 
from  0  to  90  degrees.  Figure  7 1  corresponds  to  this  case.  We  begin  with  the  breathing  wave  and 
extensional  wave  response  for  the  specially  orthotropic  structure  shown  in  figure  71(a).  As  wrap 
angle  increases,  we  observe  a  slowing  of  the  extensional  wave  phase  velocity.  The  breathing  wave 
phase  velocity  first  decreases  and  then  shows  an  increase,  with  the  change  occurring  at 
approximately  45  degrees.  This  change  in  phase  velocities  corresponds  to  a  shifting  of  stiffness 
from  the  longitudinal  to  the  circumferential  directions.  Between  the  angles  of  0  and  90  degrees, 
we  see  another  branch  developing  and  then  almost  disappearing  at  the  40-degree  position,  as  its 
phase  velocity  at  first  increases  to  approximately  the  45-degree  point  and  then  decreases  after  that. 

This  additional  branch  corresponds  to  the  v  component  of  the  n  =  0  mode  of  propagation, 
which  only  exists  when  the  shell  becomes  generally  orthotropic.  Figure  37  illustrates  the  torsional 
motion  of  this  mode,  which  occurs  because  there  is  coupling  between  the  shear  strain  and  the 
normal  stresses  (equation  (7)). 
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Figure  67.  Bending  Shell  Transfer  Surface  (n  =  0);  dB  =  lOlog ((Pj(al2)/Px)2) 


86 


rad/m  k  rad/m  k  rad/m 


Bending  Shell  Sequence  Continued 
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Figure  68.  Bending  Shell  Transfer  Surface  (n  =  1);  dB  =  10log((P,-(a/2)/Px)2) 


Effects  of  Wrap  Angle  <))  on  the  n-0  Mode 

S.  Orthotropic  n  = 


Frequency  Hz 

(a) 

G.  Orthotropic  n  = 

=  2 

Frequency  Hz 

(b) 

G.  Orthotropic  n  = 

=  4 

Frequency  Hz 

(c) 

Figure  71.  Bending  Shell  TYansfer  Surface  (n  =  0)  for  0  <  <|)  <  90;  dB  =  101og((/>i-(a/2)//>^.)'2) 
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Effects  of  Wrap  Angle  <j>  on  the  n-0  Mode  Continued 

G.  Orthotropic  n=  0  6=  6 
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Figure  71.  Bending  Shell  Transfer  Surface  (n  =  0)  for  0  <  0  <  90;  dB  =  101og((P,<a/2)/Px)2) 

(Cont’d) 
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SHELL  WALL  STRAIN  TRANSFER  SURFACES 

Thus  far,  we  have  considered  responses  in  the  inner  fluid  created  principally  by  the 
displacement  W.  Let  us  now  turn  our  attention  to  the  shell  and  explore  the  strain  occurring  in  the 
shell  wall  due  to  an  applied  radial  pressure  excitation  from  the  outer  fluid.  Additional  theoretical 
development  is  required  to  obtain  equations  for  the  strain  in  the  shell  from  the  displacement 
magnitudes,  which  are  provided  in  equation  (89). 

Theoretical  Development 

From  equation  (89),  the  displacement  magnitudes  will  be  calculated  from  a  P(,  excitation  only. 
The  stress  resultants,  equation  (13),  will  be  calculated  using  the  displacement  magnitudes.  The 
actual  stresses  will  be  calculated  from  the  stress  resultants.  The  strains  will  be  calculated  from  the 
stresses  and  then  rotated  to  the  desired  angle  H1.  It  should  be  noted  that  4*  does  not  have  to  equal  <|> 
but  can  be  an  independent  angle  of  choice. 

From  equation  (89),  the  displacement  magnitudes  of  a  point  in  the  shell  are 


and 


(SjiS-)-,  S2|6p) 

W  =  -  -----  -  P, 

det\~J 


_  (S  12^23  ^13^22)  p 

U  —  -  -  ----  r  o, 

det\s\ 


y  _  (  ^11^23  +  ^  13^21  ^  p 

det  s 


(96) 


(97) 


(98) 


The  stress  resultants  needed  to  compose  the  stresses  are  Nxx,  N qq,  /Vlfl,  Mxv,  Mee,  and  Mx$.  We  will 
use  the  following  linear  relationship  for  the  actual  distribution  of  stress  across  the  thickness  of  »he 
shell: 


G  = 

V.V 


Nxx  »2  Mxxz 


°ee 


/V0e  12  MqQ: 


°V0  = 


/VA.0  \2MxqZ 


(99) 

(100) 

(101) 
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It  should  be  noted  that 


°'e  ~  h 


"a  12M,e- 


°0v  "  h 


NQx  12  M0vr 


11  h  (102) 
will  not  necessarily  yield  identical  results.  There  is,  then,  some  objection  to  the  assumed  linear 
stress  distribution.  However,  as  Fluggels  points  out,  the  difference  is  usually  small  in  thin  shells 
and  can  be  ignored.  We  should  add  that  a  better  approach  might  be  to  take  an  average  of  the  two 
shear  stress  equations  in  equation  (102)  rather  than  use  equation  ( 101 )  alone. 

After  the  necessary  derivatives  of  the  displacements  have  been  taken  and  inserted  into 
equation  (13),  the  stress  resultants  are  written  in  terms  of  the  displacement  magnitudes  as 

N„  =  i[U(Dnk*Dt  ”)+V(D n"+Dl<tk  +  KKi)}+W{D'2+Kl 3"+*,/'}, 

cT  a  u  a~ 


n  n 


=  i{U(Di2k  +  D2ft  +/C26  ,)+l'(D22'+D26(t)} 

U  n  Q 


+  “+  --  (1  -/,-) -/^26' 3  }, 

U  St  fl 


=  nuwlbK  +  Dj'a)  +r<D*"+DMt  +  /r(i6*  ))  +wiD*  +  k,^  +  km"a}. 


=  - '  ( UK, ,  *  +  l- (K ,,  ",  +  2K * )  )  +  -IV  i  Kuk2  +  If,,";  +  2*,/* } . 

a  u  a~  a 


A -,6  ( 1  -  //“)  ->  .it 

^ee  =  /  ;  { Un  -  Vka }  +  W{K22  ,  -  Kx2k~  -  2  K2"k  } , 


a**  a 


(103) 
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Equations  (99)  through  (101)  are  used  to  solve  for  the  stresses.  The  final  steps  now  involve  the 
calculation  of  the  strain.  The  strain-stress  relations  are  then19 


[e  1 

XX 

eee 

= 

e.v0 

Sm,,  Sm,2  Sm,6 

12  Sm22  S m26 

16  Sm26  S™33 


09 


xQ 


(104) 


The  Sm  and  Sm  coefficients  are  given  by 

Sm 


u 


E  ’ 


5m  j2  =  - 


12 


Sm12  =  £  , 

^m66  =  ♦ 
u12 

Sm,,  =  Sm,,cos4<])+  (2Sm,2  +  Sm^)  sin2<(>cos2<|>  +  Sm22sin4(j>, 

Sm,2  =  Sm,2(  sin4  <))  +  cos4<j>)  +  (Smu  +  Sm22~  Sm^)  sin2(()cos2<j), 

Sm22  =  Sm,,  sin4  (j>  +  (2Sm,2  +  Sm^)  sin2<|>cos2<t)  +  Sm22 cos4 <j), 

Sm,6  =  (2Sm,,  -2Sm, 2  -Sm^)  sin<|)cos3<f)  +  - (2Sm22-2Sm, 2  -Sm^)  sin3 <pcos<t>, 

Sm26  =  (2Sm,,  -2Sm,2-Sm66)  sin3 <j)cos<|)  +  -  (2Sm22  -  2Sm]2  -  Sm^)  sin<|>cos3<t>, 

5m66  =  2(2Sm,, +2Sm22- 4Sm, 2 -Sm^)  sin2(|)cos2<t)  +  Sm66(sin4(|)  + cos4 ()>)  • 

The  strains  in  the  shell  coordinate  system  are  then  rotated  to  the  desired  angular  position,  \y,  using 
the  rotation  matrix  given  by  equation  (106): 

cos2\y  sin2\|/  -2sin\j/cos\y 
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The  primes  on  the  subscripts  of  the  stresses  refer  to  the  coordinate  system  specified  by  the  angle 
y,  where  is  measured  as  0  is  (see  figure  3).  If  =  <{>,  then  the  following  relation  exists: 


err 

£.l 

£2'2' 

II 

e22 

£1'2 ' 

_£>2 

(107) 


Transfer  Surface  £\\IPa 

There  is  a  large  volume  of  data  available  to  the  analyst  at  this  point.  The  wrap  angle  of  the 
reinforcement  <j>  can  be  varied.  The  strains  are  available  at  any  angle,  \|/,  and  at  any  location  across 
the  shell  thickness  from  The  strain,  £n,  is  the  change  in  length  per  unit  length  of 

composite  shell  material  when  the  stress  is  applied  to  the  1  face  and  in  the  1  direction.  The  1 
direction  corresponds  to  a  principal  material  property  direction. 

We  will  consider  the  following  restricted  case  consisting  of  £\\/P0 ,  calculated  at  z  =  h/2,  with 
the  simplification  that  \}/  will  be  chosen  to  be  equal  to  $.  Therefore,  we  will  be  examining  strains 
in  the  shell  along  the  <|>  direction. 

Liquid-Filled  Shell 

Let  us  begin  by  returning  to  the  urethane  shell  with  inner  and  outer  fluids.  In  figure  72,  the 
transfer  surface  for  the  longitudinal  strain  has  been  calculated.  The  two  branches  correspond  to 
the  breathing  and  extensional  waves  in  an  isotropic  shell.  Over  this  limited  range  of  wavenumber, 
the  slope  of  both  branches  is  quite  linear;  at  this  point,  the  breathing  wave  has  not  become 
dispersive. 


Frequency  Hz 

Figure  72.  Bending  Shell  Transfer  Surface  (Liquid-Filled);  dB  =  101og((£jj(/i/2)/Fo)2) 
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Air-Filled  Shell 


For  the  remainder  of  this  discussion,  we  will  change  the  properties  of  the  inner  fluid  to  those 
corresponding  to  air,  results  are  presented  as  the  stiffness  of  the  shell  and  the  reinforcing  fiber  are 
gradually  increased  for  mode  n  =  0.  The  properties  of  the  inner  fluid  are  now 

p.  =  1.2  kg/m3, 

Cj  =  343  m/s  . 

We  will  begin  with  the  transfer  surface  for  the  isotropic  cylinder  in  figure  73.  To  provide  an 
understanding  of  what  type  of  wave  propagation  each  branch  represents,  we  plot  the  individual 
displacement  magnitudes  in  figure  74  for  k  =  9. 15  rad/m. 

For  the  first  resonant  peak  in  figure  74,  U  exhibits  the  maximum  displacement.  The  second 
and  third  resonant  peaks  are  modes  of  propagation  in  w.  We  see  coupling  between  the 
displacements  as  we  did  before  in  figure  46.  There  is  also  a  strong  similarity  to  figure  66,  the  case 
of  the  soft  compliant  core.  The  properties  of  the  soft  core  are  very  similar  to  those  of  air. 
Therefore  the  internal  stiffness  of  the  system  of  figure  66  is  governed  by  the  core,  and  the  effect  of 
the  inner  fluid  has  been  diminished. 

In  the  circumferential  direction,  £22 '!?<>  *s  plotted  in  figure  75.  Figure  76  is  a  comparison  of  the 
longitudinal  and  circumferential  strains  at  k  =  0.1  rad/m.  From  this  figure,  we  see,  as  expected, 
that  the  circumferential  strain  is  some  20  dB  larger  than  the  longitudinal  strain. 

When  longitudinal  reinforcement  is  added,  the  branches  reorient,  as  shown  in  figure  77. 
Again,  we  turn  to  analyzing  the  displacements  corresponding  to  this  specially  orthotropic 
structure.  The  first  resonance  in  figure  78  now  corresponds  to  a  mode  of  propagation  in  w, 
followed  by  a  mode  in  u  for  the  second  resonance.  The  third  resonance  behaves  as  in  figure  74,  a 
mode  in  w.  Increasing  the  shell  stiffness  in  the  longitudinal  direction  has  increased  the  phase 
velocity  of  the  u  component  of  mode  n  =  0  (figure  77).  The  order  of  the  branches,  from  left  to 
right,  now  corresponds  to  w,  u,  w  components  of  mode  n  =  0,  respectively. 

For  the  generally  orthotropic  shell  (figures  79  and  80),  we  find  the  transfer  surfaces  to  include 
a  new  branch,  a  v  component  of  mode  n  =  0.  The  second,  third,  and  fourth  branches  are  now 
modes  in  which  U  is  the  maximum  displacement  (figure  80).  Figures  74,  78,  and  80  provide  a 
very  clear  reminder  that,  for  the  n  =  0  axisymmetric  mode  of  propagation,  v  =  0  for  isotropic  and 
specially  orthotropic  shells.  An  n  =  0  vibration  only  produces  a  v  displacement  in  generally 
orthotropic  shells. 
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Figure  73.  Isotropic  Bending  Shell  Transfer  Surface  (Air-Filled); 
dB  =  I01og((en(/i/2)/Po)2) 


Figure  74.  U  and  W  Component  Amplitudes  at  k  =  9.15  rad/m  From  Figure  73 
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Strain  In  The  Shell  Continued 
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Figure  79.  G.  Orthotropic  Bending  Shell  Transfer  Surface  (Air-Filled); 

dB  =  101og((en(/f/2)//>o)2) 


Figure  80.  U,  V,  and  W  Component  Amplitudes  at  k  =  15.5905  rad/m  From  Figure  79 
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Air-Filled  Shell  of  Increasing  Stiffness 


All  the  data  presented  thus  far  have  been  for  the  urethane  shell.  The  modulus  of  elasticity  for 
the  shell  matrix  and  reinforcement  will  be  increased  to  the  following  values  in  figures  81  through 
89.  The  shell  will  remain  air  backed  while  it  is  immersed  in  the  outer  fluid  (whose  properties  are 
given  with  the  data  on  the  urethane  shell). 

Shell  Data  for  Figure  81 

Em  =  1.5  x  10s  Pa,  Cm  =  0.01, 
v  =  0.40, 

Ej-  =  1.5  x  1010  Pa,  C/=0-0, 

V/  ~  0.15. 


Figures  82,  83,  and  85  Through  89 

Em  =  2.5  x  109  Pa,  £m  =  0.01, 

v*  =  0-40^ 

Ej-  =  7.0  x  1010  Pa,  £  =  0.0, 
vf  =  0.15. 

Shell  Data  for  Figure  84 

=  2.5  x  109  Pa,  Cm  =  0.l, 

=  0.40, 

Ef  =  7.0  x  1010  Pa,  ^=0.0, 
vf  =  0.15. 

Comparing  figures  79,  81,  and  82,  we  see  all  the  branches  becoming  less  dispersive  as  the 
stiffness  of  the  shell  is  increased.  The  phase  velocity  increases  as  well  with  this  progression  in 
stiffness. 

Figures  83  and  84  compare  two  levels  of  damping.  Broadening  of  the  width  of  the  branches  is 
obvious  in  figure  84.  The  resulting  decrease  in  the  peak  level  of  the  branches  is  most  apparent  in 
figure  85,  where  we  compare  a  cut  in  wavenumber  through  both  transfer  surfaces. 
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Strain  In  The  Shell  Continued 
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Figure  81.  Transfer  Surface  (Em  =  1.5  X  109  Pa,  Ef=  1.5  X  1010  Pa); 

dB  =  101og((eu(/i/2)/Po)2) 
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Figure  82.  Transfer  Surface  (Em  =  2.5  X  10  Pa,  Ef  =  7.0  X  10  u  Pa); 
dB  =  10Iog((e„(/i/2)/Po)2) 
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Extensional  Waves 


Returning  to  figure  83,  we  observe  that  the  branches  shown  for  this  n  =  0  mode  are  the  least 
dispersive  of  all  the  cases  presented.  This  shell  is  the  best  example  of  a  fulfillment  of  the 
assumptions  made  when  the  simple  expression  for  the  phase  velocity  of  the  extensional  wave 
given  by  equation  (75)  is  derived.  The  effects  of  the  fluids  and  bending  stiffness  were  neglected  in 
the  derivation  of  equation  (75)  along  with  the  frequency-dependent  terms.  The  second  branch  in 
figure  83  corresponds  to  the  u  component  of  mode  n  =  0.  From  equation  (75),  a  phase  velocity  of 
2105  m/s  is  calculated.  From  figure  83,  the  reciprocal  of  the  slope  of  the  branch  yields  2136  m/s. 
Both  values  are  in  very  close  agreement  and  provide  a  good  example  of  the  case  where  the  simple 
expression  is  valid. 

Expanding  the  frequency  range  to  1 1,500  Hz,  we  will  look  at  transfer  surfaces  for  the  stiffest 
shell  (figure  82)  up  to  mode  n  =  3  as  before.  Once  again,  the  strain  plotted  in  the  surface 
corresponds  to  the  <j)  direction  shown  above  each  surface. 

As  frequency  increases,  the  n  =  0  modes  become  dispersive.  Because  of  the  extreme  stiffness 
of  this  shell,  only  the  first  two  components,  w>  and  m,  of  the  n  =  1  mode  of  propagation  are  visible 
within  the  frequency  range  chosen  for  display.  Only  the  w  component  of  the  n  =  2  and  n  =  3 
modes  appears  in  figures  88  and  89.  The  cutoff  frequencies  are  extremely  high  for  such  stiff 
shells. 

General  observations  concerning  modes  n  =  0  through  n  =  3  reveal  an  increase  in  cutoff 
frequency  with  increasing  mode  number  as  well  as  decreasing  phase  and  group  velocity  for  a 
given  component  of  a  mode  at  a  fixed  wavenumber  for  a  given  category  of  shell  stiffness  (i.e., 
isotropic,  specially  orthotropic,  or  generally  orthotropic).  This  means  that  energy  propagates  more 
slowly  in  the  branches  corresponding  to  a  particular  component  of  a  mode  as  the  mode  number  is 
increased  for  a  fixed  wavenumber. 
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Magnitude:  dB  =  1  O*log(l£n/P  I2)  @  761  Hz 
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Figure  »5.  Comparison  of  a  Cut  Through  Figures  83  and  84  at  761  Hz 
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Strain  In  The  Shell  Continued 
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Figure  86.  Transfer  Surface  (n  =  0);  dB  =  101og((Eu(6/2)/Po)2) 
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Figure  87.  Transfer  Surface  (n  =  1);  dB  =  10log((E||(/i/2)/l>o)2) 
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Strain  In  The  Shell  Continued 
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Figure  88.  Transfer  Surface  (n  =  2);  dB  =  101og((Ejj(/i/2)/Fo)2) 
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CONCLUSIONS 


From  first  principals,  we  have  formulated  the  equations  of  motion  for  a  generally  orthotropic 
cylindrical  shell  excited  by  a  plane  wave  imparting  either  a  normal,  longitudinal,  or 
circumferential  stress  to  the  surface  of  the  shell.  The  shell  itself  can  be  prestressed  in  the  axial 
direction  by  static  tension.  The  shell  is  immersed  in  an  outer  fluid  and  contains  an  inner  fluid;  the 
two  fluids  need  not  be  similar.  The  partial  differential  equations  of  motion  are  solved  for  the  case 
of  harmonic  excitation  in  the  x-  and  0-directions.  Harmonic  solutions  are  considered  for  both  a 
standing  wave  and  a  traveling  wave  in  0.  Either  is  appropriate  for  the  simpler  cases  of  isotropic 
and  specially  orthotropic  composite  shell  material.  Only  the  traveling  wave  solution  in  0  is 
appropriate  for  the  generally  orthotropic  shell.  The  resulting  wave  propagation  is  confined  by  the 
finite  but  continuous  circumference  of  the  shell  to  an  integer  number  of  wavelengths  around  the 
circumference,  where  the  number  of  wavelengths  is  given  by  mode  number  n. 

The  shell  exhibits  coupling  between  the  three  displacements  w,  u,  and  v.  For  the  isotropic  and 
specially  orthotropic  shell,  a  radial  excitation  gives  rise  to  coupled  u  displacement  for  the  n  =  0 
mode  of  propagation,  while  for  n  >  1  both  u  and  v  coupled  displacements  occur.  In  the  generally 
orthotropic  shell,  we  find  coupled  u  and  v  displacements  for  all  modes  (/i  =  0  included).  The 
dynamic  behavior  of  the  isotropic  and  specially  orthotropic  shell  is  characterized  by  two  branches 
for  mode  n  =  0  and  three  branches  for  modes  n  >  1 .  In  the  generally  orthotropic  shell,  there  are 
three  branches  for  any  mode  number. 

Considering  the  magnitude  of  the  shell  response  as  a  function  of  wavenumber  and  frequency 
for  a  fixed  mode  number,  we  see  a  number  of  branches  exhibiting  cutoff  frequencies  for  modes 
n  >  1 .  The  cutoff  frequency  marks  the  onset  of  propagation  for  a  given  mode.  At  the  cutoff 
frequency,  the  group  velocity  approaches  zero  and  the  phase  velocity  approaches  infinity.  Energy 
propagates  at  the  group  velocity;  therefore,  at  cutoff,  no  energy  propagates  in  that  particular 
branch.  Cutoff  frequencies  occur  at  zero  wavenumber  for  the  isotropic  and  specially  orthotropic 
shell;  for  the  generally  orthotropic  shell,  there  is  an  offset  from  zero  wavenumber  at  the  cutoff 
frequency.  For  a  given  component  of  a  mode  of  propagation,  the  cutoff  frequency  increases  with 
increasing  mode  number  n.  Bending  stiffness  also  influences  the  cutoff  frequency.  The  w 
component  of  modes  2  and  3  does  not  exhibit  a  cutoff  frequency  for  the  membrane  case.  The 
bending  shell,  however,  does  exhibit  a  cutoff  frequency  for  this  mode  of  propagation. 

The  transfer  surfaces  for  the  generally  orthotropic  shell  are  asymmetric  with  respect  to  the 
wavenumber  axis.  This  shell  presents  a  different  stiffness  to  the  excitation,  depending  on  whether 
the  excitation  travels  in  the  +jt-  or  -jc-direction.  This  difference  in  stiffness  is  responsible  for  the 
asymmetry  that  results  in  the  wavenumber-frequency  surface. 

Imaginary  wavenumbers  used  in  the  excitation  and  the  solution  for  the  shell  response  result  in 
inhomogeneous  vibrations.  The  amplitude  decreases  with  increasing  x,  while  it  remains  harmonic 
and  homogeneous  with  respect  to  0. 

The  vibration  of  the  shell  creates  pressure  fields  in  both  fluids.  For  a  given  frequency,  when 
the  wavenumber  of  the  vibration  of  the  shell  surface,  k,  is  less  than  or  equal  to  the  propagation 
wavenumber  in  the  medium,  the  pressure  field  propagates  in  the  fluids.  When  the  wavenumber 
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exceeds  the  propagation  wavenumber  in  the  medium,  an  evanescent  pressure  field  results  whose 
decay  from  the  surface  of  the  shell  is  functionally  dependent  on  /„  in  the  inner  fluid  and  Kn  in  the 
outer  fluid.  The  effect  of  the  outer  fluid  on  the  structure  is  to  decrease  the  cutoff  frequency,  as  was 
observed  for  the  w  component  of  the  n  =  2  mode.  An  internal  core  slows  down  the  phase  velocity 
for  the  h’  component  of  the  n  =  0  mode  of  propagation. 

A  normal  pressure,  Pa,  excites  a  high  level  breathing  wave  and  a  weak  extensional  wave.  The 
shear  stress  excitation,  Pv,  excites  a  high  level  extensional  wave  and  a  moderate  level  breathing 
wave.  The  extensional  wave  is  the  most  nondispersive  mode  of  propagation,  as  it  engages  the 
shell  in  a  minimum  amount  of  bending  and  very  slight  interaction  with  the  fluids.  For  the  w 
component  of  the  n  =  0  mode  of  propagation,  as  wrap  angle  increases,  we  observe  a  slowing  of 
the  extensional  wave  phase  velocity  and  at  first  a  decrease  and  then  an  increase  in  the  phase 
velocity  of  the  breathing  wave,  with  the  shift  in  speed  occurring  for  a  wrap  angle  of 
approximately  45  degrees.  This  characteristic  corresponds  to  stiffness  shifting  from  the 
longitudinal  to  circumferential  directions. 

When  considering  shell  wall  strain,  the  first  step  is  to  calculate  the  stress  resultants  from  the 
displacement  magnitudes.  From  the  stress  resultants,  the  stress  distribution  across  the  shell 
thickness  is  obtained  based  on  the  linear  superposition  of  the  normal  and  bending  stress.  The  air- 
backed  shell  provides  a  good  match  for  the  simple  linearized  equation  for  the  extensional  wave 
speed  (equation  (75)).  General  observations  concerning  modes  n  =  0  through  n  =  3  reveal  an 
increase  in  cutoff  frequency  with  increasing  mode  number  as  well  as  decreasing  phase  and  group 
velocity  for  a  given  component  of  a  mode  at  a  fixed  wavenumber  for  a  given  category  of  shell 
stiffness  (i.e.,  isotropic,  specially  orthotropic,  or  generally  orthotropic).  This  means  that  eneigy 
propagates  more  slowly  in  the  branches  corresponding  to  a  particular  component  of  a  mode  as  the 
mode  number  is  increased  for  a  fixed  wavenumber. 

At  the  time  of  this  writing,  the  computation  speed  of  workstations  allows  for  the  evaluation  of 
the  stresses  and  strains  as  a  function  of  two  independent  variables.  We  have  displayed  the  results 
of  this  study  as  grey-scale  images  or  wire  frame  surfaces.  Each  image  is  composed  of  either  1 28 
by  128  points  (figure  47)  or  256  by  256  points  (figure  64).  Currently,  it  is  possible  to  exercise  the 
solutions  as  a  function  of  three  independent  variables,  resulting  in  a  volume  rather  than  a  surface 
representation.  Such  calculations  are  performed  more  expediently  on  supercomputers,  although  it 
is  becoming  possible  to  work  with  volume  visualization  on  workstations. 
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APPENDIX  A 


STRESS  RESULTANT  INTEGRATION  Nm 

The  following  integrals  are  necessary  for  the  integration  of  Nqq: 

c  x 

=  x  -  alog  (x  +  a) 


(A-l) 
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Only  the  first  two  terms  will  be  retained  in  the  series  expansion  for  the  logarithm  ratio,  which 
results  in 
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The  above  simplification  means  we  are  neglecting  terms  of  order  —  and  higher.  Substituting 

a 

equation  (10)  into  (11)  and  then  into  the  equation  for  Nqq  in  equation  (1),  and  using  the  following 
notation  for  the  derivatives,  yields 


dx 


()=()' 


)  =  (  )' 


(A-4) 


A-l 


resulting  in 
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Evaluating  at  the  limits  h/2  and  -h/2,  the  z  terms  cancel  out,  which  gives 
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The  log  terms  in  equation  (A-7)  are  simplified  according  to  equation  (A-3): 
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Collecting  terms  and  simplifying  equation  (A-8),  we  arrive  at  the  following  expression  for  AW 
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If  the  stiffnesses  of  equation  (12)  are  applied  along  with  a  reversal  of  the  shorthand  notation  in 
equation  (A-4)  to  the  above  expression  for  AW  we  will  obtain  the  expression  given  in  equation 
(13)  for  AW 
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APPENDIX  B 


DETERMINATION  OF  THE  COMPOSITE  MATERIAL  PROPERTIES 

Ej,  E2,  Vj2j  AND  Gj2 

The  determination  of  the  composite  material  properties  can  be  obtained  by  either 
experimental  measurement  or  with  an  analytic  model  if  both  the  properties  of  the  matrix  and  the 
fiber  are  known.  Presented  here  are  two  models  for  the  composite  properties:  one  based  on  the 
strength  of  materials  approach  (i.e.,  the  rule  of  mixtures)  and  the  other,  on  an  elasticity  approach.1 

Rule  of  Mixtures 

The  most  prominent  assumption  made  for  the  mechanical  behavior  of  the  composite  material 
is  that  strains  in  the  fiber  direction  (the  1  direction)  are  the  same  in  the  fiber  as  in  the  matrix 
material.  The  volume  fraction  of  fiber  and  matrix  will  be  defined  by  the  following  equation: 


The  composite  modulus  in  the  1  direction  can  be  expressed  as 


Ex  =  EfVf+EmVm 


(B-l) 


(B-2) 


In  the  transverse  direction,  the  same  transverse  stress  is  assumed  to  be  applied  to  the  fiber  and 
the  matrix,  resulting  in  the  following  for  E2: 


E  -S&L 

2  V„Ef+ 


tf+VfZn 


(B-3) 


The  major  Poisson  ratio  v12  is  given  by 


V 12  =  Vn,y'm+Vf';r 


(B-4) 


The  corresponding  Poisson  ratio,  v21,  is  found  from  the  condition  of  symmetry  of  the  compliances 
to  be 


v 


21 


(B-5) 


1.  R.  M.  Jones,  Mechanics  of  Composite  Materials,  Hemisphere  Publishing  Corporation, 
New  York,  1975,  pp.  90-114. 
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The  assumption  for  the  shear  modulus,  Gl2,  is  that  the  shear  stresses  on  the  fiber  and  the 
matrix  are  the  same,  resulting  in 


G12  - 


V/ 


ymGf+Vp, 


(B-6) 


Elasticity 

One  very  useful  elasticity  solution  to  the  composite  problem  is  presented  by  Tsai2  for  the  case 
in  which  the  composite  is  represented  by  a  combination  of  fibers  in  contact  with  each  other  and  by 
fibers  interspersed  uniformly  in  contact  with  the  matrix.  The  two  states  are  depicted  in  figure  B-l . 
Most  composites  are  a  combination  of  both  states  and  would  be  represented  by  a  linear 
combination  of  the  two  solutions  with  the  value  of  C  ranging  from  0  to  1 . 


2.  S.  W.  Tsai,  Structural  Behavior  of  Composite  Materials,  National  Aeronautics  and  Space  Administration.  CR-71, 
July  1964. 
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Tsai  obtains  the  following  elasticity  solution  to  the  problem  when  contiguity  is  considered: 


E2  =  2(1 -v/+  (vf-vm)VJ 


(1  ~C) 


+c 


(Kf(2Km  +  Gm)-Gm{KrKm)Vm) 
VKm  +  Gm)+2(Kf-Km)Vm 

(Kf(2Km  +  Gf)  +  Gf(Km  -  Kf)  Vm) 
VKm  +  Gf)-2(Km-Kf)Vm 


with 


*/  = 


2(1-V/)’ 


Km  = 

m 


G(  = 


2(l-vm)’  /  2(l+v/)’ 


m  ”  2(1+  v  )  ’ 

m 
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(B-8) 


12 


=  (1-0 


+  c 


Kpf(2Km  +  GJ  V/  +  Knym  (2Kf+  GJ  Vn 
Kf(2Km  +  Gm)-Gm(KrKm)Vm 

(2 Kf+  Gf)  Vm  +  Kfvf(2Km  +  G,)  V, 


KfV*m  +  Gm)+Gf{Km-Kf)Vm 


and 


G 


12 


,  2Gr  (GrGJvm  .  rr  (G/+GJ  -  «r-Gm)vm 

-20,+  (Gf-  Gm)  Vm  '(Gf+Gm)  +  (GrGJVm 
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The  modulus  in  the  direction  of  the  fibers  is  modified  from  the  rule  of  mixtures  to  account  for 
fiber  misalignment.  The  misalignment  factor,  ka,  which  varies  from  0.9  to  1 .0  and  is  dependent  on 
the  manufacturing  process: 

E\  =  M^/+W  •  m-in 
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APPENDIX  C 


EFFECTS  OF  AN  ELASTIC  SOLID  CORE  ON  THE  INNER 
PRESSURE  FIELD,  FOR  THE  AXISYMMETRIC  CASE 


Let  us  consider  the  interaction  between  the  pressure  field  in  the  inner  fluid  and  an  elastic  solid 
core.  The  core  will  be  centered  on  the  axis  of  the  cylinder  as  depicted  in  figure  C-l. 


Figure  C-l.  Elastic  Core  Cross  Section 


The  core  radius  is  rc  with  Young’s  modulus  Ec  and  Poisson  ratio  v(..  The  core  is  modeled  as  a  long 
thin  cylinder  using  the  plain  strain  approximation.  The  mass  of  the  core  will  be  neglected  in  this 
development  and,  as  stated,  the  results  will  be  limited  to  the  axisymmetric  case.  The  plane  strain 
simplification  is  most  appropriate  here  and  will  be  used  to  describe  the  elastic  behavior.  Equation 
(C-l)  shows  the  variation  of  the  stresses  with  radius,  and  equation  (C-2)  defines  the  strain-stress 
and  strain  displacement  relations:1 


C 


c  =  2C,  + 

'  2  2 

r 


3 


(C-l) 


and 


1.  C.  Wang.  Applied  Elasticity ,  McGraw-Hill  Book  Company.  New  York,  1953.  p.  55. 
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1  +  V  1  +  V 

er=  -g— ^  1(1  -vr)or-vco9],  Ee=  -g-1  l-vcor  +  ( 1  -  vc)  <Je] , 


duc 

e=-r, 

r  dr 


e„  =  “i 


(C-2) 


For  a  solid  cylinder,  C3  must  be  zero  for  the  stresses  to  be  finite  at  r  =  0.  Because  C3  is  equal 
to  zero,  we  find  that  ar  =  ae  =  ic2  =  c.  Together  with  equation  (C-2),  we  have  the  following 
equation  for  the  displacement  in  the  r  direction: 


Uc  = 


(1  +vc)  (1-2 \c)rC 


(C-3) 


The  velocity  potential  in  the  inner  fluid  will  now  contain  both  Bessel  function  terms. 
Referring  back  to  equation  (39),  we  have 

<P,  =  {AJ0(lr)  +BY0(lr)}ei(kx-<0,) , 

,2  0)2  ,2  f  ,2  ®2 

/  =  — 5-  -  k  for  k  <  —  • 


(C-4) 


The  velocity  and  pressure  field  for  the  fluid  are  given  according  to  equations  (36)  and  (37)  and 
result  in 


w,.  =  {-  M7,  (Ir)  -lBYl(lr)}eiikx-<0'),  Pj  =  i(0Pi{AJo(lr )  +BY0  (/r)  } 


(C-5) 


The  constant  C  is  obtained  by  considering  the  boundary  condition  for  radial  stress  at  the 
surface  of  the  core,  or  =  -p  at  r  =  rc.  Because  pressure  is  a  compressive  load  applied  to  the  core, 
its  sign  is  negative,  resulting  in  the  following  equation: 


C  =  -i<nPi{AJ0(lrc)+BY0(lrc)}e 


i(kx-U)t) 


uc=  I- 


(l+vf)  ( 1  - 2vf ) r 


jz'G)P/-  {AJa  (/rc)  +  BY0(lrc)}  e 


i(kx-(Ot) 
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The  velocity  in  the  core  is 


u,  =  -0) 


p  r^{AJ0{lrc)+BY0(lrc)}e 


i  (kx  -  to/) 


(l+vr)(l-2vr) 

K 
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The  condition  of  contact  at  the  fluid-solid  interface  requires  that  the  velocity  of  the  core  be 
equal  to  the  velocity  of  the  fluid  at  r  =  rc.  Equating  equations  (C-5)  and  (C-7)  yields 


2  V, 


2„  ..  Vc, 


+Ul(lrc)}  +B{-0>‘p.rc^Yo(lrc)  +  /F,(/r(.)}  =  0 


(C-8) 


Another  equation  containing  the  undetermined  constants  A  and  B  results  from  equating  the 
fluid  and  shell  velocity  at  r  =  a(.  For  this  axisymmetric  case,  the  displacement  of  the  shell  in  the 
radial  direction  is  given  by  equation  (58);  the  magnitude  of  the  shell  velocity  is  then  «  =  -m>w, 
resulting  in 


-iwW  =  -AIJ]  (la;)  -BIYX  ( la, ) 


(C-9) 


Configuring  equations  (C-8)  and  (C-9)  into  matrix  form,  we  have 


IJ]  (la;) 


IY ,  (la;) 


v  i  r  v  i  A  =  '“H/ 

-<D2p,.r,.^y„  «r,.)  +  /■/,(/,•)  W  L° 


-  V 

..  <•, 


(C-10) 


Solving  equation  (C-10)  for  A  and  B  and  then  substituting  into  equation  (C-5)  gives  the  following 
expression  for  the  internal  pressure  field: 


to  2  2  05  2 

for  the  range  —=  >  k  with  /  =  —r  -  k  , 

ci 


Pi(r)  = 


P/(0 2{CY0(lr) -DJJlr)  }  Wei{kx-w) 

lJET^)-CY^laJ} 


Pi  (at)  = 


Pj(°2  {CYa  (la;) -  DJJla ;)  }Wei(kx-'*,) 
l{DJ\(la;)-CYi(lai)} 


P;W  { CY()  (Ir)  -DJ  n  (Ir)  } 


p.o r{CY()(la;)  -DJJIa,)  } 


a.(r)  =  ----- = - = - - ,  a. (a.)  =  - — =— -  - - 

'  /  {DJ ,  (Id;)  -  CY ,  (Id;)  }  '  ''  I  {DJ ,  (Id;)  ~  CY ,  (/fly)  } 


co  p  r  .V 

C  =  UlUrc)-—j—Lj0(lrc). 


_  CO  p  /  •  V 

D  =  /F,(/r(.)--^— ry„(/r.)  . 
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In  a  similar  fashion,  for  the  corresponding  region  of  wavenumber 


k 2 


with 


Pi(r) 


P/.Q)2  {Hl0  ( mr )  -  GK0  ( mr )  }  ^ 
m{GKx  (ma{)  +///,  (ma,)  } 


p,.co2  {///0  (ma,)  -  GK0  (mat)  >  HV(far_aM) 

p,(a.)  =  - - - = - , 

'  '  m{GKx  (mat)  +HIi  (««*)  } 

Pi®2  {  W0  (wr)  -  GK0  (mr)  }  p,©2  {HIQ  (ma,)  -  GKa  ( ma ,)  } 

m~{GK{  (ma^  +HI{  (majj  }  ’  ”  m  {Gtf,  (ma,)  +  W,  (ma,)  }  ’ 


rc©2p (mrc)  -  mAT,  (mrc) 
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The  effect  of  the  core  is  not  easily  understood  from  equation  (C-12).  Small  aigument 
approximations  to  the  Bessel  functions  will  be  used  to  simplify  the  expression  for  a,  in  order  to 
observe  the  effect  of  the  core  on  the  breathing  wave  speed  within  certain  limiting  approximations. 
The  small  argument  approximations  are 


/„(*)-  1,  Ka(x)  — In  (x)  ~  2, 

lx  (x)  ~  ^ ,  Kx  (x)  -  I  • 

^  X 


(C-13) 


Using  equation  (C-12),  we  arrive  at  the  following  simplification  for  a ,/a,)  for  the  case  where 


.  ®. 
k  »  —  . 
c. 


,2  r  ^  i  i,2. 


« /(«,-)  = 


-p  ,©M-+rrc} 


k(r2-a2)  v  (l+*2ah 


(C-14) 


Considering  the  case  for  a  rigid  core  ^  =  o,  a/a,),  equation  (C-14)  reduces  to 


(C-15) 
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Following  the  same  development  for  the  breathing  wave  speed  that  resulted  in  equation  (70),  we 
have  the  expression  below  for  the  case  of  a  rigid  core: 


The  specially  orthotropic  stiffness  matrix  in  equation  (4)  simplifies  equation  (C-16)  to  the 
following  expression,  while  making  use  of  the  symmetry  property  of  equation  (73): 


It  is  clearly  seen  that  the  presence  of  the  core  slows  the  breathing  wave  down  relative  to  the  result 
obtained  in  equation  (74). 
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